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ABSTRACT 

] We propose and evaluate the feasibihty of a new strategy to search for planets 

0^ I via microlensing observations. This new strategy is designed to detect planets 

in "wide" orbits, i.e., with orbital separation, a, greater than ~ 1.5Re- Planets 
>■ ' in wide orbits may provide the dominant channel for the discovery of planets 

o ■ 



> 



via microlensing, particularly low-mass (e.g.. Earth-mass) planets. Because the 
ongoing microlensing observations and extensions of them should be able to 
discover planets in wide orbits, we provide a foundation for the search through 
detailed calculations and simulations that quantify the expected results and 
(^n [ compare the relative benefits of various search strategies. If planetary systems 

similar to our own or to some of the known extra-solar systems are common, 
then the predicted detection rates of wide-orbit events are high, generally in the 
range 2 — 10% of the present detection rate for apparently single events by stars. 
^ I The expected high rates should allow the microlensing observing teams to either 

place significant limits on the presence of planetary systems in the Galactic 
Bulge, or begin to probe the population in detail within the next few years. We 



2 ' also address the issues of (1) whether planets discovered via microlensing are 



O ■ 

^ ' likely to harbor life, (2) the feasibility of follow-up observations to learn more 

about planet microlenses, and (3) the contamination due to stellar populations 
^ of any microlensing signal due to low-mass MACHOs. 

c3 ' Subject headings: - Gravitational lensing: microlensing, dark matter - Stars: 

planetary systems, luminosity function, mass function - Planets & satellites: 
general - Galaxy: halo - Methods: observational - Galaxies: Local Group. 

1. Introduction 

When planets are discovered through monitoring programs similar to those presently 
being carried out, the distance to the lens will typically be on the order of kiloparsecs. 
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Although some spectral follow-up may be possible (§6), imaging, such as that planned 
to probe possible companions to nearby stars (Stahl & Sandler 1995; Bender & Stebbins 
1996; Labeyrie 1996; Angel & Woolf 1997) will not be possible in the near future. 
Nevertheless, planets discovered as microlenses can play an important role in developing 
our understanding of planetary systems in our Galaxy and beyond. The reason for this 
is precisely because microlensing observations probe vast volumes of the Milky Way and 
other galaxies, such as the Magellanic Clouds and M31. Thus, microlensing provides a 
unique window for studying the statistics of planetary systems (numbers and properties) 
and their dependence on the local stellar environment. Microlensing searches complement 
velocity-based searches in another way as well. Particularly for the wide systems studied 
here for the first time, the microlensing searches can be effective for low-mass planets 
orbiting at low speeds. 

The framework for the work to date on discovering planets through microlensing was 
established by Mao & Paczyhski (1991) and Gould & Loeb (1992). These authors found 
that when the separation, a, between the star and planet is close to the Einstein radius. Re, 
of the star (O.Si?^; ^ a ^ 1.5-R^;), a significant fraction of events (~ 5 — 20%) in which the 
star serves as a lens would be perturbed in a detectable way. This has been referred to as 
"resonant" lensing, both because the separation must be close to -R^; in order for the signal 
to be detectable, and because the signal itself is sharp and distinctive. Q The detection 
strategy is to monitor light from an ongoing event at frequent intervals in order to observe a 
short-lived perturbation. With the analysis of more than 200 microlensing events reported 
to date, it is not clear whether any planetary-lens events have been discovered. Indeed, 
the fraction of microlensing events to date in which evidence of a planet in a resonant 



■^We will use the phrase "zone for resonant lensing" to refer to the region in which distinctive events 
occur, whose characteristics tend to be related to the presence of the caustic structure (i.e., events of the 
type specifically studied by Gould & Loeb [1992], Bennett & Rhie [1996], and Wambsganss [1997]). The 
spatial extent of the zone for resonant lensing has been studied most recently by Wambsganss (1997); our 
results do not rely heavily on the exact position of its boundaries as a function of mass ratio. 

^ One candidate for a resonant planetary lensing event has been suggested (Bennett et al. 1996). There 
are, however, two reasons to be cautious about the interpretation of this event. First, the mass ratio derived 
from the binary-lens fit is ~ 0.043, which is large enough to be consistent with lensing by a binary stellar 
system. Second, the degeneracy of the physical solution has not yet been worked out. The degeneracy 
may be of two types: (a) Other binary solutions with values of q differing from the one for this fit by as 
much as an order of magnitude may prove to be equally good fits (Di Stefano & Perna 1996); this needs 
to be systematically checked, (b) Other physical effects may prove to be important, such as finite source 
size effects (Witt & Mao 1994; Witt 1995), finite lens-size effects (Bromley 1996), and blending (Di Stefano 
& Esin 1995, Kamionkowski 1995). (In fact blending has been a feature of every other binary-lens event; 
Udalski et al. 1994, Alard, Mao, and Guibert 1995, Alcock et al. 1997a.) 
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orbit is detectable is apparently less than ~ 1%. This is not necessarily inconsistent with 
the original predictions, however, because planets with the appropriate distance from the 
central star may not be common. In addition, finite-source-size effects can wash out the 
signature of the short-time-scale perturbations characteristic of planetary mass ratios, 
decreasing the number of detectable events (Bennett & Rhie 1996). 

The purpose of this paper is to develop a complementary framework for the discovery 
of planets by microlensing. The basic idea is that planets located more than ~ 1.5Re 
from the central star can give rise to microlensing events that reveal evidence of their 
presence. We will say that planetary orbits in which the separation between the star and 
planet is larger than 1.5 Re are "wide". Concrete definitions as well as our expectations 
for the possible existence of a population of wide-planetary-orbit lenses are discussed in §2. 
Because the ongoing observational programs have already developed the tools needed to 
begin a systematic search for planets in wide orbits, the rest of the paper is designed to give 
a detailed view of what observers can expect to see and to compare different observational 
strategies. In §3 we focus on the character of the events expected when the lens is a wide 
planetary system. In §4, we present the results of Monte Carlo simulations of lensing by 
several model planetary systems, including our own solar system, to test basic features of 
detectability. In §5 we study the link between the search for evidence of planets in wide 
orbits and the search for extraterrestrial life. Blending may play an important role in 
allowing us to learn more about planetary systems, including the mass of a planet and 
the spectral type of the central star; blending is the topic of §6. When the size of the 
source is comparable to or, to some extent, larger than the Einstein ring of the planet, the 
probability of detecting planets in wide orbits can increase. This is in contrast to lensing 
by planets located in the zone for resonant lensing, where there is a tendency for detection 
probabilities to decrease as the ratio between the source radius and the Einstein radius 
increases. Q Thus, planets in wide orbits may provide the primary channel for the discovery 
of low- mass (e.g.. Earth-mass) planets. Finite-source-size effects are the subject of §7. In 
addition to including multiple planets, planetary systems may be more complex, containing 
moons bound to planets, and also belts or clouds of asteroids and comets. Because the 
detectability of these features tends to be limited by finite-source-size effects, we discuss 
them in §8. §9 is devoted to a discussion of the detection strategies that will improve our 
ability to discover planets in wide orbits. In §10 we address the likely near- and long-term 



^ A more precise statement is that for a planet of a given mass, there is a range of source sizes for which 
finite-source-size effects influence the rate of detection, but do not prevent the detection of planet-lens events. 
Within this range, finite-source-size effects tend to increase the rate of detection of planet-lens events when 
the planet is in a wide orbit, and to decrease the detection rate for an identical planet located in the zone 
for resonant lensing. The net effect is to increase the relative rate of detectable wide-orbit lensing events. 
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results of implementing those strategies. Finally, in §11 we summarize our results. 

2. Wide Planetary Systems 
2.1. Definition 

We will say that a planet is in a "wide" orbit if its distance from the central star is 
large enough that the isomagnification contour associated with A = 1.34 is comprised of 
two separated, closed curves, one centered near the star and the other centered near the 
planet. For mass ratios ranging from 10^^ to 10^'^, the critical orbital radius, aw, beyond 
which an orbit becomes wide is roughly equal to 1.5Re- It is also useful to define an inner 
radius, Og. For a < ac the fraction of all light curves displaying deviations from the standard 
point-lens form is small (generally less than 1%). For Oc < a < a larger fraction of light 
curves exhibit perturbations, including caustic crossings. The value of Uc is approximately 
equal to O.SRe for values of q appropriate to planets (Gould & Loeb 1992). 

In Figure 1 we consider planetary lenses in close, resonant, and wide orbits; in each case 
the mass ratio between the planet and the star is 0.001, comparable to that between Jupiter 
and the Sun. Thousands of events were generated in which the track of the source passed 
behind each lens; in this way we could sample the results expected when the system serves 
as a lens. In the plots of distributions of event durations, the most striking feature is the 
appearance of a peak corresponding to short-duration events when the orbital separation is 
wide (right-most panel). Indeed, 3% of the events for a = 3Re have durations less than 
2 days. Lensing by stellar systems produces such a distinct peak of short-duration events 
only when there are wide-orbit planetary systems in which planets serve as isolated lenses. 
Other features illustrated by the figure include the following. When a planet is in a close 
orbit, the probability of detecting its presence via microlensing is small. At "resonant" 
separations, the fraction of events revealing evidence of the caustic structure associated 
with the presence of the planet is significant. While the fraction of such events decreases as 
the orbital separation increases toward the "wide" regime, a significant fraction of lensing 
events in the wide regime are (1) the short-duration events mentioned above, which are 
absent in the other cases, and/or (2) multiple-peak events, in which the separation between 
two of the peaks is so large, and the magnification between them falls so low, that they are 
best thought of as repeating events. 



- 5 - 




Fig. 1. — Plotted is the distribution of event durations for 3 microlenses. Each lens has a 
planet of mass mj in orbit with a star of mass M; in each case we chose q — rrii/M to he 
10~^, the mass ratio between Jupiter and the Sun. The lenses differed from each other only 
in the value of the orbital separation, a. For each lens, several thousand source tracks were 
randomly generated, and the light curve was computed for each event, fee represents the 
fraction of all events which exhibited caustic crossings; fmuiu if 2, /s) represents the fraction 
of all events which exhibited multiple (2 ,3) peaks. 
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2.2. Extending Planetary Searches Beyond the Resonant Case 

The excitement about resonant lensing by planets was fueled in part by a wonderful 
coincidence. Gould & Loeb (1992) noted that, if a system identical to our solar system 
happened to be located halfway between our position and the center of the Bulge, and if the 
system were viewed face-on, the separation of the planet corresponding to Jupiter from the 
system's star would be very close to the value oi Re associated with the mass of this star. 
That is, Jupiter would be in the zone for resonant lensing. There are, however, two reasons 
to be cautious about using this example to limit the search for planets to those in resonant 
orbits. First, we do not know that the spatial relationship between the Sun and Jupiter is 
an example of a universal property of planetary systems. Mindful of this, Bennett & Rhie 
(1996) have constructed a "power-of-2" planetary model, in which the separation between 
each planet and the central star increases by a factor of two for each successive planet. In 
such a model, most planetary systems can be expected to contain one planet in a resonant 
orbit. 

A second reason to avoid limiting the microlensing searches to resonant planets is 
that the value of the Einstein radius depends not only on the stellar mass, but also on the 
relative positions of both source and lens to the observer. 

\AMGDs x{x-l)Y 

where M is the mass of the lens, Ds is the distance from the observer to the lensed source, 
and X = Dl/ Dsi with representing the distance from the observer to the lens. 

In Figure 2 we consider planets with orbital separation a = 2aw = "i Re- For these 
systems we show the relationship between the stellar mass, M, and a;, when the orbital 
period is fixed. We have assumed that the orbital plane is the same as the lens plane, and 
show those values of x for which both lens and lensed source are located in the source 
galaxy (the Bulge, the Magellanic Clouds, or M31). If the stellar mass is in the range from 
0.1 to IOM0, the orbital periods of planets in wide orbits range from a few years to a few 
hundred years, with larger values more typical for more distant galaxies. 

In addition, the projected value of the separation depends on the angle, a, between 
the normal to the orbital plane and our line of sight. Consider an orbit with a > a^. Let a 
represent the angle between the normal to the plane of the orbit and the normal to the lens 
plane. For a < amin, where 

ttmrn = cos~^(a^/a), (2) 

the projected value of the orbital separation is always greater than a^- If o; > amin, then 
for a circular orbit, the fraction of time for which the projected separation is greater than 
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Fig. 2. — Each curve corresponds to a fixed value of the orbital period; shown are curves 
for 1 yr, 10 yrs, 100 yrs, 1000 yrs. Each curve illustrates the relationship between the stellar 
mass, M, and x for a planet in a wide orbit, with a = 3 Re- In each panel the range of the 
variable x corresponds to lenses located within the same galaxy (Bulge, Magellanic Clouds, 
M31 [left to right]) as the lensed stars. Note that the range of masses and orbital periods 
best explored via wide lensing events is different for different source galaxies. 
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Qyj (so that the planet can be viewed as wide) is 



/ = - cos 
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(3) 



If the orbit is eccentric, then the separation between the planet and star will be wide for an 
even greater fraction of the time. 

Thus, even if planetary systems had uniform properties, the fact that they are located 
at different spatial positions and are tilted at different angles relative to our line of sight 
would still not favor resonant orbits over others. Bennett & Rhie's power-of-2 model, or 
something like it, is necessary to ensure that most planetary systems have one planet in the 
zone ttc < a < Qw associated with resonant lensing. This is because a^/ac ~ 2; thus, if the 
i*'* planet has separation Oj from the central star, and the power-of-2 system is inclined so 
that Oj cos a ~ Oc, we have Oj+i cos a ~ a^. 

The known planetary and brown-dwarf binary systems are considered in Table 1, 
together with a "power of 2" model and a "power of 3" model. The known systems included 
in the table were selected from the Encyclopedia of Extrasolar Planets (accessible from 
www.obspm.fr/darc/planets/encycl.html). Two criteria were used: (1) if the planetary 
system were placed in the Bulge and viewed face-on, the planet (or brown dwarf) listed 
had to be in an orbit that would either be wide or located in the zone for resonant lensing, 
(2) the existence of the planet or brown dwarf needed to be hsted as "confirmed". In the 
power-of-2 and power-of-3 model, the first planet listed is the innermost planet that would 
be in a wide orbit, were the planetary system to be placed in the Bulge and viewed face-on. 
We note that the next planet inward would have a good chance of being viewed in the zone 
for resonant lensing. Indeed for both of these theoretical models, there is some ambiguity 
associated with the arbitrarily-chosen position of the closest planet to the central star. 
Because of this, we carried out a set of calculations in which we averaged over the position 
of the innermost planet. We found that, if we truncate the radius of the planetary system 
at ~ 10^^ cm, there are on average ~ 9 — 10 planets in wide orbits for every planet located 
in the zone for resonant lensing. To derive the numbers shown in the table, we averaged 
over inclination angle a, and used equations (2) and (3) to determine the fraction of time 
any given planet would be in a wide orbit. As above, we have assumed that the orbits are 
circular. 

The considerations in this section show that planets in wide orbits exist, and may 
be common. In fact the number of planets in wide orbits may be as much as an order of 
magnitude larger than the number of planets in resonant orbits. The relative contribution 
of these two classes of planet lenses to the detection rate depends on our abihty to identify 
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the associated events. 



3. The Detection of Wide Planets 

There are two distinctive signatures of the presence of wide planets in the stellar 
systems that serve as lenses. The first is a significant number of short-duration events in the 
data sets. Although relatively few individual short-time-scale events may be unambiguously 
identifiable as being due to planets, the statistics associated with the presence or absence 
of such events provide a strong diagnostic for properties of the population of planetary 
systems. The second distinctive signature is a small number of repeating events. Repeating 
events are good diagnostics for individual planetary systems, yielding a value for the mass 
ratio and the orbital separation (expressed in units of Re) as projected onto the plane of 
the sky. 

For both isolated and repeating events, if the central star is luminous and provides a 
significant fraction of the light received along the line of sight to the lensed source, then 
spectroscopic studies, conducted during the event and also at baseline, can determine the 
central star's spectral type and may allow us to derive the mass of the planet (§6). Another 
diagnostic of the planet-lens mass for both isolated and repeating events is provided by 
finite-source-size effects, which can, in addition to providing information about the mass of 
the planet-lens, also help to increase the detection rate of low-mass planets in wide orbits 
(§7). 



3.1. Notation: Encounters, Events, and Repetitions 

The planetary systems we consider consist of a central star and N planets moving 
in wide, circular orbits. Let the index i label the central star (i = 0) and the planets 
{i — l,N). The characteristics of the planetary system are specified by the mass of the 
central star, M = mo- the planetary masses, rrii, the associated mass ratios, Qi — rrii/M, 
and the orbital radii, a^, of the planets. For the sake of convenience, we will sometimes drop 
the index for the star; unless otherwise specified, M and Re refer to the mass and Einstein 
radius of the star. 



-10- 



Table 1. Planets in wide and resonant orbits in known and model systems; 

Ds = 10 kpc, X = 0.9. 



Planetary system Planet a (i) %dose '^^^ %res %Me. '"^^ 



Solar system: 


Jupiter 


5.2 


5.9 


21.0 


73.1 




Saturn 


9.5 


1.7 


4.6 


93.7 




Uranus 


19.2 


0.4 


1.0 


98.6 




Neptune 


30.1 


0.2 


0.4 


99.4 




Pluto 


39.8 


0.1 


0.2 


99.7 


Known extrasolar systems: 


55 Cnc 


4.0 


11.0 


89.0 


0.0 




HD 29587 


2.5 


37.0 


63.0 


0.0 




PSR B1620-26 


38 


0.1 


0.3 


99.6 




Gl 229b 


40 


0.1 


0.2 


99.7 


Theoretical power-of-2 system: 


1 


4.8 


7.0 


27.0 


66.0 




2 


9.6 


1.6 


4.5 


93.9 




3 


19.2 


0.4 


1.0 


98.6 


Theoretical power-of-3 system: 


1 


4.8 


7.0 


27.0 


66.0 




2 


14.4 


0.7 


1.9 


97.4 




3 


43.2 


0.1 


0.2 


99.7 


•^^^ Separation of the planet from the central star in astronomical units (AU). We consider circular orbits. See the text for the 



criteria used to select the planets listed here. 

(■^'Percentage of the time during which the planet is "close" (a < 0.8/?^;), averaged over all inclinations. 

(3)Perccntagc of the time during which the planet is "resonant" (0.8-Rs < a < 1.5-Re), averaged over all inclinations. 

Percentage of the time during which the planet is "wide" (a > 1.5-Rb), averaged over all inclinations. In the power-of-2 
and powcr-of-3 models, all planets beyond the third had a probability greater than 99% of being located more than 1.5 Re from 
the central star. 



- 11 - 



3.1.1. Width of the Lensing Region 

Let Amin be the minimum value of the peak magnification for which a perturbation is 
reUably ascribed to microlensing. Each value of A„jj„ is associated with a specific value of n, 
the number of Einstein radii equal to the minimum necessary distance of closest approach. 
If Amin is 1.58 (1.34, 1.06, 1.02), then n = 0.76 (1, 2, 3). Note that n need not be an integer. 

For any given lens, the "width" of the lensing region is defined to be 

Wi^nRE,i. (4) 

Thus, if Amin = 1-34, the width associated with each mass is simply equal to its Einstein 
radius. With this definition of the width, the probability of an encounter with a particular 
lens can be taken to be 2^1 = 2ni?E_j. 

3.1.2. Encounters and Events 

When a planetary system with several planets in wide orbits serves as a lens, the track 
of the source may pass through the region of infiuence of several lenses. When the source 
track passes within Wi of lens i, we will say that an "encounter" is underway. We will use 
the word "event" to refer to a source track (and the associated light curve) that experiences 
one or more encounters. 

3.1.3. Repeating, Isolated, and Overlap Events 

When more than one encounter occurs, we will dub the event a "repeating" event; 
events with just a single encounter will be called "isolated". 

The relatively short time duration of encounters in which a planet serves as a lens 
means that frequent monitoring is needed to detect them. If a stellar-lens encounter takes 
100 days, a Jupiter-mass lens encounter will last ~ 3.3 days, and an Earth-lens encounter 
will last ~ 8 hours. (Finite-source-size effects can, however, increase the duration of 
detectable encounters. We consider them in §7.) Thus, without frequent monitoring, we 
may miss the occurrence of an isolated planet-lens event. If one component of a repeating 
event involves the central star, then we may see that component, and miss the planet-lens 
component. To increase our chances of discovering repeating events, we can use the fact 
that the encounter due to the star is likely to occur first roughly half the time, and continue 
frequent monitoring of lensing light curves, even after they have returned to baseline. The 
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calculations we describe in §4 indicate that ~ 50 — 75% of the cases in which a planet 
will subsequently serve as as a lens can be discovered if frequent monitoring continues 
for 100 days. Similarly, moons revolving about planets could be discovered with frequent 
monitoring lasting for a time interval on the order of days, subsequent to a short-duration 
encounter which might be due to lensing by a planet (§8). 

The ideal case is that in which the perturbation due to lensing by one mass has not yet 
ceased before lensing by the second mass begins to have a noticeable effect. We will refer to 
such events as "overlap" events. Overlap events can occur when the separation, a, between 
the lenses is less than wi + W2. As a decreases below wi + W2, many events which involve 
encounters with both lenses will give rise to light curves in which the magnification does not 
fall below baseline between encounters. These light curves are similar to others exhibiting 
repetitions, except that the perturbation due to the first encounter is still noticeable when 
the second encounter begins. We are almost guaranteed to discover 1/2 of all repeating 
"overlap" events, even with no or little change in detection strategy. As the value of a 
decreases further, a larger fraction of the encounters with the smaller lens are part of events 
in which the larger lens is also encountered; for a < Wi — W2, this is true of all encounters 
involving the smaller lens. Most of the associated light curves will exhibit a single connected 
perturbation, with the structure in the wings largely determined by the most massive lens, 
and the structure near the peak largely determined by the less massive lens. Although 
such light curves may appear to be highly anomalous, we should be able to fit them to the 
appropriate lensing model. Since the light curve perturbations are long-lived (relative to 
encounters with an isolated planet), we might be able to detect and identify 100% of such 
overlap events. (See Figure 3.) 

In general we expect at most -|- 1 encounters and N repetitions. In principle, if the 
source travels slowly and the planets orbit rapidly, a single planet could give rise to more 
than one encounter. In practice this is extremely rare (§3.4.6). Repetitions, however, are 
not rare. As we will see, if most stars have planetary systems similar to our own, a data set 
containing the number of events already observed along the direction to the Bulge should 
contain one or more events in which a repetition is due to lensing by a wide orbit planetary 
system. As we will also discuss, the detection of such events can be optimized by frequent 
sampling with sensitive photometry. 

3.2. Normalization of Event Rates 

We will be interested in the rate at which detectable events occur, and how that rate 
is influenced by changes in the detection strategy. The Einstein radius of the central star 
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= source track 




2.0 



Fig. 3. — The innermost planet of the power-of-2 system. The dotted circles represent the 
1.0-, 1.5-, and 2.0-Re rings for the central star. The large solid circle is the orbit (viewed 
face-on) of the innermost wide planet, at a separation of 4.8 AU. If the central star is of 
solar mass, then the radius of the small solid circle on the orbit is roughly equal to the 
Einstein radius of a 10- Jupiter-mass object. When wq — Re — wo,o, there is a range of 
angles along which a source track can pass through the Einstein ring of the planet without 
also encountering the star. When wq — 2 Re — 2^0,0, any encounter with the planet also 
results in an encounter with the star. 
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provides a convenient normalization. The rate of encounters in which the central star serves 
as a lens, with the magnification due to the central star achieving a peak value of at least 
1.34 {Amin = 1-34), is proportional to 2 Re- Because the lensing width of the central star 
(when Arain = 1-34) will play an important role, we define wo,o to be equal to Re- If the 
photometric sensitivity is such that Amin can be smaller than 1.34 (i.e., n > 1), so that 
the effective width of events involving the central star is larger (wq > Wo,o), the rate of 
all detected events, including those due to the central star, will increase. Nevertheless, 
2 ^0,0 is a convenient normalization constant, because it allows us to compare the rates we 
compute to the presently-measured rate of events. This is as follows. Along the direction 
to the Bulge, events are being discovered by the MACHO team at a rate of roughly 50 
per year (See §9 for the distribution of observed durations.) The present detection criteria 
used by the MACHO team are strict: A^in = 1.58, corresponding to n (in Eq. 4) being set 
equal to ~ 0.76. In addition, other teams, surveying different fields with somewhat different 
strategies, are also surveying the Bulge. For example, the OCLE team, which in its first 
incarnation discovered ~ 12 Bulge events, has recently brought a new telescope on-line 
(Udalski, Kubiak & Szymanski 1997; Paczynski et al. 1997). It is therefore reasonable to 
assume 75 — 100 events of the type we use for our normalization per year along the direction 
to the Bulge. Thus, because we use 2 wq.o to normalize our results, when we find that a 
particular detection strategy leads to a rate of detectable events of a certain type (e.g., 
events with 2 repetitions) equal to p%, this means that between 0.75 p and p such events 
could be discovered per year along the direction to the Bulge. 



When a wide planetary system serves as a lens, the most common type of event is one 
in which the track of the source passes through the Einstein ring of a single mass. Because 
the rate of encounters associated with a given lens scales as the square root of its mass, 
encounters in which the star serves as a lens are the most frequent. Using the normalization 
described above, the rate of lensing due to planets is approximately given by 



Because isolated events are generally 1 — 2 orders of magnitude more common than repeating 
events. Pi generally provides a good estimate of the rate of isolated events involving planets. 
In the next subsection, we will focus on the rates of repeating events, which should be 
subtracted from the above expression to make it exactly correspond to the rate of isolated 
events. The large majority of short-duration isolated events will be indistinguishable from 
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Isolated Events 




(5) 



1=1 
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the standard point-mass-lens light curve. A small number will exhibit caustic crossings. 

The time duration of isolated planetary-lens events, relative to the time duration of 
stellar-lens events, also scales as ^/g^. The implication is that a distribution of events due 
to lensing by stars with wide-orbit planets is necessarily accompanied by a distribution of 
shorter-duration events. The fraction of events in the latter distribution is proportional to 
the average value of ^g, and the position of the peak or peaks also provides a measure of 
the mass ratios typical of planetary systems. 

3.3.1. Identifying Peaks In The Distribution of Short Events 

A peak corresponding to short-duration encounters in the distribution of encounter 
durations should be unambiguously identifiable as being due to the presence of lenses 
with small masses. We argue as follows that such a peak is unlikely to be hidden by the 
short- duration tail of encounters due to lensing by stars. Three circumstances can make a 
stellar-lens encounter have a short measured duration. 

(1) The relative velocity Vt an observer would measure between the source and the lens 
may be exceptionally large. Since, however, the distribution of transverse velocities for 
planet-lens and stellar- lens encounters should be roughly the same (see Eq. [T6| ), we should 
be able to disentangle velocity effects, at least for statistical samples of encounters. 

(2) The track of the source may just graze the lens' Einstein ring. In this case, the peak 
value of the magnification allows us to measure 6, the distance of closest approach [which 
we will also call Dmin], and to determine that, although short, the observed encounter was 
nevertheless due to the presence of a lens with a larger Einstein ring than expected for a 
planet. These corrections are implemented by the observing teams. 

(3) Light from the lensed source could be blended with light from other sources along 
the line of sight. If, however, the effects of blending are dramatic enough to shorten the 
duration of an event significantly, then we should be able to detect evidence for the blending 
and to subtract its effects from the event before comparing its duration with that of other 
events. Because blending can provide a valuable tool for the study of planet-lenses, we 
discuss it in more detail in 56. 



3.4. Repeating Events 
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3.4- i- One Repetition 

When the source track passes close to one wide planet, and also within the region of 
influence of the central star and/or that of other wide planets, the event will appear to 
repeat. Such repeating events were studied by DiStefano & Mao (1996) in the context of 
wide binaries-i.e., binaries in which the separation between the components is large enough 
that the A = 1.34 isomagnification contour consists of two disjoint closed curves. 

Let Wi {W2) represent the width associated with lensing by the more (less) massive lens. 
We assume that the separation between the two lenses is wide-i.e., the isomagnification 
contours associated with A = 1.34 are distinct. As we discussed in the previous section, this 
generally means that if one lens is the central star and the other is a planet, a > 1.5 Re- f\ 
The widths, wi and W2, however can be smaller or larger than the orbital separation, since 
their values are tied to issues of photometric sensitivity. 

With 1^0,0 defined as above, the detection rate for repetitions involving lenses 1 and 2 
is given by: 



2 1 f 

^1,2 = - 7^ { Omax{Wl + W2) + 6min{w2 - Wi) + a 

Vr 2 Won I 



C0s{9max) " COs{ 



(6) 



In this expression, the values of 6max and 6min depend on the value of the orbital separation, 
a, as compared to the widths, wi and W2. 



For a > wi + W2, 



For Wi + W2 > a > wi — W2, 



sin ^{—^—^);9min = sin \— —). (7) 



"max 1 "min ollt y j. yd j 

For a <Wi — W2, 

9=9=- (9) 

"max "mm V"^/ 

The limit 'Wi + W2«a corresponds to the case in which the separations are extremely 
wide. In this case, 



n,2 = -'^ = n^ 

TT a Won 



2 Re,1 Re,2 
TT a Won 



(10) 



^Dynamical stability of the planetary system implies that two planets will generally be separated from 
each other by significantly more than 1.5 times the Einstein radius of the more massive planet. 
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The rate of repetitions is inversely proportional to the separation. Another key feature of 
this expression is the quadratic dependence of the width of the lensing region, which is 
directly related to the photometric sensitivity. 

The second extreme limit corresponds to the case when the separation between the 
lenses is smaller than Wi. In this case, 



^1, 



W2 



n 



R 



E,2 



Wo,0 



That is, the rate of repeats involving lenses 1 and 2 is the same as the rate of events 
involving 2; all source tracks passing through the lensing region associated with 2 necessarily 
pass through the lensing region associated with 1. 



3.4.2. Estimates 

For the purposes of microlensing, the most important difference between a stellar 
binary and a planetary system is that the planetary system may contain several planets. A 
naive generalization of Eq. 10 leads to the following expression for the average probability 
of repeating events. 



2n^ 

irTsys Jo 



N 
i=l 



di{t)\ 



N-l 

+ EE 

1=1 j>i 



di{t)-aj{t)\ 



dt. 



(12) 



where Tgys is the time taken for the configuration of the planetary system to approximately 
repeat, and di(t) is the position vector of the ith planet, expressed in units of the stellar 
Einstein radius. Note that we are considering the regime a » Wi + W2; thus the rate of 
overlap events is not included in these calculations. 

To develop a feeling for the numbers, we will consider a simple model, a power- 
of-A; model, in which the plane of the planetary system is coincident with the lens 
plane and each planet has the same mass ratio q with the central star. Using 
Jf)"^" dt\dj{t) — di{t)\^^ ~ (Oj + 0^)"-*^/^, allows us to simplfy equation ( [T^ ) in a way that 
does not overestimate the rate of wide-orbit lensing events. 



vrai 



N-l 

E 

i=0 



N~2 

+ y v- , 



(13) 



where ai is the distance between the star and the first wide planet. For the power-of-2 
model the value of the first term in the above equation is approximately 2, and that of the 
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second term is roughly 1.86y/g. The rate of repeating events in which the star and one 
planet each serve as a lens is (4n^yg)/(7r ai) . If g = 0.001 and ai = 2, this becomes 
0.02 72^. [] The observed ratio between the rate of repeats involving both the star and one 
planet, to the rate of stellar-lens encounters is 0.02 n. The ratio of such repeats to isolated 
planet-lens encounters is roughly (2 n)/(7r a A^) , where N is the number of planets in wide 
orbits. Note that the contribution of repeating events saturates at moderate values of N, 
while the rate of isolated short-duration events is proportional to N. The ratio of repeating 
events in which both encounters are due to planet lenses, to those in which one component 
is due to the central star, is suppressed by a factor 0{y/q). 



3.4-3. Multiple Repetitions 

Multiple repetitions can occur when the stellar system contains 3 or more masses. To 
derive analytic expressions, we first consider lensing by a static planetary system whose 
orbital plane is aligned with the lens plane. The effects associated with the planets' 
velocities and with changing the orientation of the orbital plane are then briefly considered. 
The analytic expressions allow one to make intuitive predictions. In §4 we turn to numerical 
simulations to derive detailed results that can be checked against the predictions. 



3.4-4- Multiple Repetitions in the Static, Face- On Approximation 

If a source track crosses through the orbit of a planet, the probability that it will cross 
within Wi of the planet itself is approximately [(2 Wj)/(7r aj)]. (We assume that the planet's 
motion during the time the source crosses its orbit can be ignored, and that curvature 
effects are also unimportant.) Define Px to be the probability that the source track will 
encounter a specific subset, X = {/5i,/?2, • • • ,Pk}, of the system's masses, producing an 
event with k — 1 repetitions. Considering the regime a » Wi + wj, the probability of 
encountering these objects, and only these objects, is 

px^^ n (^)- (14) 

if the central star is one of the masses encountered. If the central star is not encountered, 
then a similar expression is derived, but the contribution is suppressed by a factor on the 



^We note that, if n approaches 2, then a planet in an orbit with a = 2 would be detected through overlap 
events; a would need to be larger for the approximation a >> wi + W2 to apply. 
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order of We note that the curvature of the innermost planet in the set X may be 
significant and can increase the probabihty that a track will encounter this planet. The 
probability of an event with k — \ repetitions is obtained by summing over all distinct sets, 
X. 

The expression above is useful because they clarify the functional form of the probability 

of observing a particular combination of k objects: generally a rough proportionality to 
the ratio between the width of each planet's lensing region and its distance from the 
central star. This makes it simple to estimate the effects of varying the parameters of the 
source-lens system. 

3.J1..5. Orbital Inclination 

In general the plane of the orbit will be inclined relative to the line of sight. Let a 
represent the angle between the normal to the lens plane and the normal to the orbital 
plane. An approach analogous to the one sketched above for the face-on case can be used. 
The polar symmetry of the system is broken, however, since the projection of an inclined 
circular orbit on the sky is an ellipse; the probability of multiple objects acting as lenses 
depends on the direction from which the source approaches the lensing system. For some 
directions of approach, the range of impact parameters leading to detectable events is now 
smaller, but a larger fraction of events will involve encounters with more than one object, 
and so will appear to repeat. The effect can be intuited by multiplying each factor of in 
the above equations by the geometrical factor 



I / — 9 

V V sin 9 + cos^ 9 cos^ a 

where 6 is the angle the track of the source makes with the semimajor axis of the elliptical 
projection of any planet's orbit. Averaging over all possible angles of approach, there is a 
net increase in the number of repeating events. 

As a tends toward 71/ 2 and we view the planetary system edge-on, most directions 
from which the source can approach can lead only to isolated events. For a small swath 
of tracks, however, the probabihty that the source will pass through the Einstein radii 
of several planets approaches unity. If the planets are of equal mass and events can be 
detected when the distance of closest approach is Wi — nRE,i, then the probability of an 
event in which the central star and all of the wide orbit planets serve as lenses is roughly 
nRE,i/D, where n is given by Eq. 4, Re^i is the Einstein radius of a planet, and D is the 
distance between the two most widely separated objects in the planetary system. The 
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Table 2. Predictions: wide-orbit lensing events for known and model systems 



System (2) 


Pl-P© 


psun (4) 


pnosun (5) 


psun (6) 


pnosun (' 


Known systems: 












Gl 229 


19.9 


0.7 


0.0 


0.0 


0.0 


PSR B1620-26 


10.0 


0.4 


0.0 


0.0 


0.0 


Model systems: 












Power-of-2 


26.0 


2.2 


0.037 


0.017 


8.2 X 10-5 


Power-of-3 


14.7 


1.7 


0.013 


7.1 X 10-3 


1.2 X 10""' 


Power-of-4 


12.6 


1.5 


7.1 X 10-3 


4.5 X 10-3 


3.7 X 10-^ 



^^"^ The computations that produced the results in this tabic were based on the analytic approximations discussed in §3. Only 
encounters with wide planets were considered; the results were averaged over all inclinations of the system to the line of sight, 
and over all angles of approach for the source track. The detectability threshold is A^i„ = 1.34 (separation of I.ORe)- No 
minimum event duration was required for detectability. 

(^'55 Cnc and HD 29587 are not listed here because, as can be seen in Table 1, the orbital separations of the planets in these 
systems do not fall into the regime for wide lensing for the values of Dg and x under consideration. 

(3) Percentage of isolated events in which a planet is encountered. 

Percentage of events in which two encounters occur, and in which the central star is one of the objects encountered. 

Percentage of events in which two encounters occur, and in which neither object is the central star. 

Percentage of events in which three encounters occur, and in which the central star is one of the objects encountered. 

Percentage of events in which three encounters occur, and in which none of the three objects is the central star. 

(*'ln each of these systems, the planet closest to the sun is in an orbit of radius a = 4.8 AU. The system is cut off at a 
maximum possible separation of 10'' AU. Thus the power-of-2 system has 12 planets, the power-of-3 system has 7 planets, and 
the power-of-4 system has 6 planets. 
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probability of encounters with all but one or all but two of the masses in the planetary 
system are of roughly the same size, typically 0(0.001). Thus, when the microlensing teams 
have carefully followed the progress of thousands of events, they will have found several in 
which the source track traces the global structure of a planetary system that gives rise to 
several detectable encounters. 

The orbital inclination also influences the wait times between encounters. Averaging 
over all possible directions of approach should cause a net decrease in the wait times 
between encounters for repeating events, and should also produce a dispersion in the 
distribution of wait times. 

3.4.6. Velocity Effects 

The planets that can be detected via microlensing are generally far enough from the 
central star that their orbital speed, Vorb, is low. If a = /iRe then 

Vorb = 7 - 

s 

This shows that the planets that can be discovered through microlensing, and especially 
wide planets, tend to be orbiting with fairly low speed. Hence, the transverse speed Vt of 
the source (relative to the central star of the lensing system) is likely to be large enough 
compared to Vorbi that the static approximation used above can provide a good guide to 
the event probabilities. Nevertheless, there will be a small number of events in which the 
relative magnitudes and orientations of the transverse source velocity and a planet's orbital 
velocity will influence the encounter probability or the characteristics of an observed event. 
This can occur when Vt is drawn from the low-velocity end of the velocity distribution, 
particulary if x is close to unity (or zero) and the central star is massive. For example, if 
the transverse velocity of the source is small enough that RE,i/vt is a significant fraction of 
the orbital period of planet i, the event probability increases. This is because the source 
spends a nontrivial amount of time within Re^i of the planet's orbit; if the planet does not 
lie in the source's path at the beginning of this interval, it may move to cross the source's 
path during the interval, resulting in a detection. In addition, event durations can be 
significantly decreased or increased, depending on the angle between vt and Vorb- 




(16) 
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4. Simulations of Lensing by a Planetary System 

To include all of the the relevant effects, we have carried out sets of simulations in which 
specific planetary systems serve as lenses. Each simulation was a numerical experiment in 
which the projections onto the lens plane of a large number of randomly selected source 
tracks were followed. The model planetary systems we have considered are (1) the solar 
system; (2) a power-of-3 model with seven Jupiter-mass planets; (3) a power-of-2 model 
with twelve Jupiter- mass planets. In all cases we neglected "inner" planets, i.e., those 
whose physical spatial separation from the central star was less than ~ 1.5 Re, where Re is 
the Einstein radius of the star. All of the orbits have been taken to be circular, with the 
motion governed by Newton's laws. For each model planetary-lens system, we have placed 
the source population in the Galactic Bulge {Ds = 10 kpc), and the lens at xDs = 9 kpc. 
The central star in each of these systems was chosen to be of solar mass. 

At the time when the source started moving along the track with some velocity, Vt, we 
started the planets in motion, each at a randomly chosen orbital phase. As time progressed, 
we tracked the position of the source and of each of the planets. We wanted to determine 
when the source track passed close enough to any mass in order for there to be a potentially 
observable "encounter" , roughly how long each such encounter would last, and how many 
encounters there would be as the source traveled along a specific track. Because the details 
of the magnification as a function of time are not needed to derive this information, we used 
only the value of the projected separation between the source and each planetary-system 
mass to determine whether an encounter was in progress. Specifically, we asked whether 
the projected separation was smaller than some pre-selected value, Wi = uRe^i- As before, 
i labels the masses in the planetary system, with i = corresponding to the star, and 
i ranging from 1 to iV for the planets. The values of n to be used were selected at the 
beginning of each simulation, [n = 1 (n = 2) corresponds to a magnification of 1.34 (1.06).] 

10^ tracks were used to sample the events expected when the solar system served as a 
lens, and 3 x 10^ tracks were used to sample the power-of-3 and power-of-2 models. Because 
even the larger numbers of tracks used for the power-of-3 and power-of-2 models do not 
provide sampling equivalent to that used for the solar system, we have normalized the 
results so that the effective linear density of source tracks was the same for each planetary 
system. We have also smoothed the power-of-2 and power-of-3 distributions to remove 
effects due to poor samphng. We note, however, that the finer features of the statistical 
distributions, i.e., those which become apparent only as the number of events increases, are 
best seen in the plots of the solar system distributions. 

We have used our simulations to derive in detail what happens for each planetary 
system, and also to test the effects of varying some of the event parameters and detectability 
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criteria. 



4.1. Event Parameters 

The characteristics of the hght curves are determined by the characteristics of the 
event. In particular, we must consider (i) the orientation angle, a, between the plane of the 
orbit and the lens plane, and (ii) the transverse speed, Vt, of the source with respect to the 
central star of the planetary-system lens. We have therefore carried out some simulations 
in which wc have varied a and Vt, in order to test their influence on the results. In more 
realistic simulations, a was chosen uniformly over the interval — 2 tt, and Vt was chosen 
from a Gaussian distribution centered at 150 km/s, with width equal to 50 km/s. 

4.2. Detectability Criteria 

As presently implemented, the search for resonant planets is carried out in two steps. 
First, the monitoring teams (e.g., EROS, MACHO, OGLE) identify microlensing event 
candidates. They monitor the flux from O(IO^) stars regularly. Some of the monitored 
fields are not visited every night; others are the targets of regular nightly monitoring and 
may occasionally be re- visited even within a single night. If the monitoring teams discover 
that an otherwise non-variable star's flux has increased significantly above baseline in 
several consecutive measurements, they typically issue an alert so that other observers can 
monitor the star more frequently. Follow-up teams have been formed to take systematic 
advantage of this opportunity (Udalski et al. 1994; Pratt et al. 1995; Albrow et al. 1996). 
Under favorable conditions, the follow-up teams carry out hourly monitoring with good 
photometry (~ 1%). 

Within this framework, there are two key elements of detectability. First, does an 
encounter last long enough to be detected? Even in principle, an encounter cannot be 
detected unless it is caught in progress during at least one observation. Reliable detection 
generally requires the event to last long enough to span at least the time interval between 
two or more consecutive observations. In some of our simulations, we have assumed that a 
minimum event duration of one day is needed in order to reliably detect the first encounter 
between a source track and the Einstein ring of a lens. This criterion can be achieved 
by the monitoring teams in some fields, using their present observing strategy. In others 
simulations, we have dropped the requirement of a minimum duration for the first event. 
This more relaxed condition is appropriate to the follow-up teams; since they can achieve 
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hourly monitoring, even an event lasting 8 hours (such as one hkely to be due to an 
Earth- mass planet) , can be readily identified if the follow-up teams attempt to discover new 
events. 

The second key element of detectability is provided by the value of the peak 
magnification: what is the minimum peak magnification, Amin, required in order to reliably 
determine that an event occurred? When the observing teams started, they tentatively 
chose Amin = 1-34, corresponding to a distance of closest approach equal to Re- It turned 
out, however, that some apparent events with magnification above 1.34 but less than 
~ 1.58 were due to stellar variability; this has led the MACHO team, for example, to use 
Amin = 1.58. It is likely, however, that this condition can be relaxed as the continued 
study of the same fields over time will allow for better identification and tracking of stellar 
variability and decrease its possible contamination of our count of true microlensing events. 
In principle, the value of Amin is set by the photometric precision of the monitoring system. 
If the photometry is good to the 1 — 2% level, then smaller values of Amin are achievable. 
Amin = 1-06 corresponds to a distance of closest approach approximately equal to 2.0Re, 
and Amin = 1-02 corresponds to a distance of closest approach approximately equal to 
3.0Re- (To achieve this latter value would require better photometry than is typical of even 
the present-day follow-up teams.) We have carried out two types of simulations. In the 
first, we have assumed that event identification was being done by the monitoring teams; 
in these, we assumed that the distance of closest approach needed to be at least as small 
as Re- In the second, we have assumed that event identification was being done by the 
follow-up teams; in these, we assumed that the distance of closest approach needed to be 
at least as small as 2 Re- We note that finite source size effects can make events detectable 
when the distance of closest approach is even larger. We return to this point in §7; the 
results presented in this section were derived under the assumption that the size of the 
lensed source could be neglected. 

We label the three sets of detectabihty criteria used in our simulations: "A" , "B" , and 

"C" . 

Criteria A: The first encounter must exhibit magnification greater than A — 1.34 
(source- lens separation less than I.ORe) for at least 1 day in order for the lensing event to 
be detected. After the detection of a first encounter, subsequent encounters can be detected 
when A > 1.06 (source-lens separation less than 2.0Re), and are not subject to minimum 
duration requirements. 

Criteria B: All lensing encounters are detected when the sourcc-lcns separation becomes 
smaller than 2.0Re- The first encounter must have a duration of at least 1 day in order 
to be detected, but after one encounter has been detected, there is no minimum duration 
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required for the detection of subsequent encounters. 

Criteria C: All lenses are detected at a source-lens separation of 2.0Re, and no minimum 
duration is required in order to detect any encounter. 

Because we want to make contact with observations, we have defined encounter 
durations and wait times to reflect the actual required monitoring times during and between 
events. The duration of an encounter is always defined to be the interval of time during 
which the source was within 2 Re^i of a lens. Wait times are defined as the time intervals 
between encounters. 

We have used the normalization of §3.2. The primary difference is that we express the 
rates in terms of percentages. Pi is the percentage of source tracks that cross through the 
Einstein radius of a single object. All percentages are computed with respect to the number 
of isolated stellar lens events (with Apeak > 1-34) in our simulations. We also frequently 
refer to these percentages as "rates", since, as discussed in §3.2, they should be roughly 
equal to the number of events of the given type expected per year along the direction to the 
Bulge. 

(Pi — Pq) is the percentage of source tracks passing through just one object, excluding 
cases in which the central star is encountered. Pi^overiap is the percentage of source tracks 
which in which the influence of encounters with 2 lenses is clearly visible in the light curve, 
but the light curve exhibits just one continuous perturbation (i.e., once the magnification 
falls below the detectability limit, it does not again rise above it). Pi is the percentage of 
source tracks crossing through the Einstein radii of i objects. Pi^overiap is the percentage of 
source tracks passing through i + 1 Einstein rings, but in which the presence of two of the 
lenses (the central star and the innermost planet in a wide orbit) is detected in an overlap 
encounter; the magnification rises above and falls below the detectability limit only i times. 



4.3. Results 
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Fig. 4. — Distribution of encounter durations for isolated events. The duration is defined 
to be the time interval during which the projected separation between the source track and 
the lens is smaller than 2 Re^z, where Re^i is the Einstein radius of the lens. The planetary 
system is a power-of-2 system with five Jupiter-mass planets; the orbital plane was coincident 
with the lens plane {a — 0). Top panel: Vt — 10 km/s. Bottom panel Vt — 80 km/s. The 
dashed curves show the duration of overlap events. 
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4.3.1. The Effects of Systematic Variation of Vt 

We used a truncated (5-planet) power-of-2 system to explore the influence of changing 
Vt on the rate of events and characteristics of events. f\ In the simulations whose results 
are shown in Figure 4, we took a = 0; i.e., the orbital plane coincided with the lens plane. 
Holding a fixed, we systematically varied Vt. The distribution of encounter durations for 
Vt = 10 and 80 km/s are shown in Figure 4. A comparison between the two cases clearly 
illustrates that, in keeping with the predictions, the overall rate of detected events is larger 
for all types of events when vt is smaller. This effect is most pronounced for repeating 
events which involve several planets, since the probability of a repeating event behaves like 
a product of detection probabilities for each planet separately. The percentage P2,overiap 
of overlap doubles demonstrates this well, since an overlap double involves at least three 
objects, two of which must be planets. P2,overiap is larger for Vt = 10 km/s by a factor of 
~ 1.8. When Vt is comparable to the orbital velocity, as it is in the top panel, there is a 
clear dispersion in the distribution of encounter durations corresponding to short- duration 
encounters in which a planet serves as a lens. This dispersion is due to the influence of 
the planets' orbital velocities; encounters can be lengthened or shortened depending on the 
angle between the planet's orbital motion and the transverse source velocity. In contrast, 
for Vt = 80 km/s, the peak due to planets is sharp and well-defined. 

As a increases, the projection of the orbits onto the lens plane become ellipses. The 
projected orbital speed along the semi-major axis is the same as before, but the projected 
orbital speed along the transverse direction is smaller. Thus, along some directions of 
approach the effects associated with the finite size of Vt become less pronounced. 

4-3.2. Systematic variation of orbital inclination 

To systematically test the results of changing a, we again used the 5-planet power-of-2 
system, this time keeping vt fixed at 80 km/s. The results for 4 simulations (a = 0°, 45°, 75°, 
and 90°) are shown in Figure 5. For a = 0°, four peaks in the distribution of wait times 



^Note that the exact normahzation described in §3.2 requires dividing by the number of encounters, 
rather than by the number of isolated events in which the central star serves as a lens. Because these two 
different ways of computing the rates lead to similar results (using the number of events yields rates that 
are ~ 2 — 10% higher), and because the number of events provides a more straightforward comparison, we 
have chosen to divide by the number of isolated stellar-lens events. 

^The reason for the truncation was simply to achieve better sampling, i.e., a higher density of source 
tracks. 
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Fig. 5. — Distributions of durations of and wait-times between encounters for events with 
one repetition. The model planetary system serving as a lens was the 5-planet power-of-2 
system. 



-29- 



between encounters are sharp and clearly visible; these correspond to wait times between 
encounters with the central star and encounters with one of the four outermost planets. 
(The innermost planet is so close that it serves as a lens only in overlap events.) The 
even spacing between the peaks is a signature of the power-of-2 model, becoming about 
0.3 ~ logio 2 f*^^ planets with >> Re- As the inclination increases, the dispersion 
mentioned in §3. appears; note that increasing inclination can only decrease wait times, 
so that each peak is "smeared" out to the left until at high inchnations a wedge-like 
shape is achieved. The detection rates for non-repeating events decreased with increasing 
inclination, but the relative detection rates for repeating events increased, as predicted. 
For a — 45°, the innermost planet's orbit comes within 0.7Re of the star, placing it in the 
zone for resonant Icnsing most of the time. The detection rates for doubles and overlap 
singles decrease dramatically because of this, but they rise again at a = 60° as the second 
planet's projected orbit becomes small enough for overlap events to occur. At a = 90° 
(edge-on), the motion of the planets brings all of them into the zone for resonant Icnsing, 
or even closer, part of the time. This decreases overall detection rates for planets in wide 
orbits. The detection rates for doubles decrease dramatically; this is because all planets he 
along the same line, and events which might have been doubles actually become triples, 
quadruples, or higher-order events. In fact, we found that for a — 90°, the percentages 
P3, P4, and P5 were all approximately 0.3%. Thus, as mentioned in §3, a relatively large 
fraction of events can exhibit multiple repetitions, as the track of the source sweeps across 
the ecliptic, crossing through the Einstein ring of several planets and the central star. 



4.3.3. The Effects of Changing the Detectahility Criteria 

To better learn how to optimize the returns from the microlensing observations, we 

tested the effects of varying the dctcctability criteria, from the most conservative set of 
criteria (set A), to the most inclusive (set C). The results are shown in Table 3, and in 
Figures 6-10. 

The planetary systems serving as lenses are our solar system, the twelve-planet 

power-of-2 system, and the 7-planet powcr-of-3 system. Figures 6 and 7 show the duration 
distributions for non-repeating events obtained when using detection criteria A and C, 
respectively. Figures 8 and 10 show the duration and wait time distributions for events 
with one repetition when using detection criteria A and C, respectively. Figure 9 shows the 
cumulative distributions for the graphs plotted in Figure 8. 

The detection criteria of set A are the most restrictive. When our solar system serves 
as a lens, with Vt = 150 km/s, only Jupiter and sometimes Saturn are able to produce 
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Table 3. Simulation Results: wide-orbit events for known and model systems 



Detect (1) 


_ Pq (2) poverlap (3) 


P2 


poverlap (5) 
2 


P3 (6) 


poverlap (7) 
3 


Solar system, 


V = 150 km/s, a = 0°: 










A 


0.3 2.2 


0.5 


0.01 


0.0 


0.0 


B 


1.8 6.0 


1.4 


0.04 


3.0 X 10-4 


0.0 


C 


4.4 6.0 


1.7 


0.05 


2.4 X 10-^ 


0.0 


Solar system, 


V = 150 km/s, a = 75°: 










A 


0.1 1.7 


0.8 


0.02 


1.8 X 10-3 


1.5 X 10-* 


B 


0.3 2.4 


1.7 


0.04 


4.2 X 10-3 


1.5 X 10-4 


C 


1.5 2.4 


2.3 


0.04 


0.01 


1.5 X 10-4 


Solar system. 


V = Gaussian, a = uniform: 










A 


0.3 2.1 


0.7 


0.02 


1.4 X 10-3 


0.0 


B 


1.4 4.2 


1.8 


0.04 


3.2 X 10-3 


0.0 


C 


3.1 4.2 


2.1 


0.05 


5.5 X 10-3 


1.5 X 10-4 


Power-of-2 system, V = Gaussian, a = uniform: 










A 


23.7 2.2 


3.0 


0.03 


0.03 


0.0 


B 


53.8 4.0 


6.6 


0.08 


0.06 


0.0 


C 


55.9 4.0 


6.6 


0.08 


0.08 


0.0 


Power-of-3 system, V = Gaussian, a = uniform: 










A 


13.0 1.9 


1.5 


0.03 


0.01 


0.0 


B 


29.3 3.6 


3.1 


0.06 


0.02 


0.0 


C 


30.4 3.6 


3.1 


0.06 


0.02 


0.0 



Descriptions of the detection criteria can be found in the text. All probabilities are given as percentages of the number of 
events in which the central star was the only lens encountered, and in which the magnification reached at least Ami„ = 1.34. 

Percentage of isolated (non-repeating) events (one "peak" in the light curve). Only a single planet-lens was encountered. 

(3) Percentage of non-repeating events which exhibited evidence of lensing by two masses. These are "overlap" events; in these 
cases the two lenses were almost always the central star and the innermost planet. 

(4) Percentage of events with one repetition (two "peaks") These are not overlap events; two well-separated masses were 
encountered. In most cases these two masses are the central star and the second planet out, but there are other contributions 
as well. 

(5) Percentage of single-repetition events in which one component consisted of overlapping encounters; the repetition was due 
to lensing by a third mass. 

Percentage of events with two repetitions (three "peaks"); all lenses were well separated, with no overlap. 

Percentage of double-repetition events in which one component showed evidence of overlapping encounters; the two 
"repetitions" were due to lensing by two other well-separated masses. If the linear density of source tracks passing through the 
power-of-2 and power-of-3 systems had been the same as for the solar system, there would have been events in these categories. 
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Single-Peak Events: v, = Gaussian, a = uniform, detect = A 
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Fig. 6. — Distributions of durations for isolated (non-repeating) events using the detection 
criteria of set A. In each panel, the peak on the left is due to encounters with planets (for the 
solar system, mainly Saturn with a few isolated events from Jupiter), and the peak on the 
right is due to the central star. Overlap events, produced almost exclusively by encounters 
involving the innermost wide planet and the central star, arc shown by the dotted lines. 
The truncation of the large-duration end of the overlap distributions in the power-of-3 and 
power-of-2 models is due solely to the poorer statistics achieved in the simulations of these 
models; all of the distributions were normalized to the linear sampling achieved for the solar 
system. 
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Fig. 7. — Distributions of durations for isolated (non-repeating) events using the detection 
criteria of set C. Here, because an encounter did not have to have a minimum duration in 
order to be detected, the peaks tend to have low-duration tails which can run together. The 
right-most peak in each panel is due to the central star of the system, with overlap events 
involving the innermost wide planet being shown by dotted lines; the peaks on the left are 
due to planetary encounters. The central feature in the top panel is a superposition of peaks 
corresponding to Uranus, Neptune, and Saturn; the encounters due to Saturn are visible as 
a shoulder at about 1.6 days. 
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events sustaining a magnification of = 1.34 for longer than 1 day. Any other planets 
would only be detected by the follow-up teams, i.e., only if a larger object (almost always 
the Sun) was encountered first. Thus, even though a realistic velocity distribution will help, 
wide-orbit planets which are slightly less massive than Saturn are difficult to detect using a 
strategy based on a set of criteria similar to that of set A. The power-of-2 and power-of-3 
models we examined contained only Jupiter-mass planets, so each planet was detectable, 
except for the small fraction of encounters in which the source track just grazes the Einstein 
ring of the planet. We note, however, that the criteria of set A are not powerful tools for 
the discovery of solar systems like our own. 

Using the detection criteria of set B (detecting even a first encounter at A = 1.06, 
instead oi A — 1.34) generally more than doubled the detection frequencies in our most 
realistic simulations, in which vt was chosen from a Gaussian distribution and a was chosen 
from a uniform distribution. Planets could be detected at a lower peak magnification, and 
more planets in the system were able to sustain a magnification of A — 1.06 for the 1-day 
minimum duration required of the first encounter. For the solar system, Saturn was often 
detected as an isolated lens, or as the first lens encountered, and Uranus and Neptune made 
occasional appearances as isolated lenses. Note that an increase by a factor of 2 is expected, 
since, even with the set of detection criteria A, n was already equal to 2 for the second 
encounter. Deviations from the factor of 2 increase are mostly associated with events which 
include encounters with a low-mass planet, and are therefore seen primarily in the solar 
system; these deviations are due to our use of a minimum duration for the detection of the 
first encounter. 

Finally, the detection criteria of set C (removing the requirement that the first 
encounter last for at least a day) made Uranus and Neptune more regularly detectable as 
isolated lenses, whereas before they were primarily detected as part of repeating events 
involving a larger mass. These planets produce significant new structure in the duration 
distribution of Figure 7. Even Pluto produces its own peak at 10~^ days (~ 10 minutes), 
although this is too short to be detected in practice. 

To assess the relative benefits of the three detection strategies used, we note that 
the largest increase in the rates of detectable planet-lens events was realized by switching 
from the set of criteria A, in which Amin for the first encounter is 1.34 to set B, in which 
Amin for both encounters is 1.06. The gains in switching from B to C, which eliminated 
the requirement that the first detectable encounter have a duration of 1-day, could be 
substantial only if, like our own solar system, the lens planetary system contains wide-orbit 
planets less massive than Jupiter. 

In general, the smaller the value of A^i^, the longer the duration of the observed 
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event anyway. Thus, while frequent monitoring is desirable, and more frequent monitoring 
than is presently achieved by the monitoring teams would be very valuable, the key issue 
is sensitivity to encounters in which the peak magnification may be smaller than 1.34. We 
note that this result is likely to remain valid, even when other effects are considered. For 
example, both blending and finite-source-size effects decrease the peak magnification. Thus, 
detection rates for events subject to these effects may also be improved if encounters which 
achieve only smaller values of can be detected. 



4.3.4. General Features 

Perhaps the most interesting feature of the results is that they clearly indicate the 
feasibility of the search for planetary systems containing planets in wide orbits. For the 
most optimistic model, the 12-planet power-of-2 model, even the strictest detection criteria 
yielded a rate of 4.5% for events showing some evidence of two or more objects in the 
system and, 23.5% for isolated short-duration events. For the detection criteria of set C, 
10% was the rate of events showing some evidence of two or more objects in the system, 
and ~ 56% was the rate of isolated short-duration events. It is clear that, if the power-of-2 
model is realized in nature with any great frequency, we will be able to observe many 
planet-lens events even during the next few years. Conversely, an absence of large numbers 
of interesting planet-lens events, particularly if we use the detection criteria of set C, would 
allow us to definitively falsify the hypothesis that most stars are accompanied by power-of-2 
planetary systems. 

While it may not be surprising that such a radical model is verifiable or falsifiable, the 
relative ease with which planetary systems such as our own can be discovered is indeed 
worthy of note. The rates of wide orbit events range from ~ 3% for the detection criteria 
of set A to ~ 10% for the detection criteria of set C. This means that, even over the 
short-term, observing strategies such as the ones we propose in §9, will lead to significant 
results. In §10 we discuss the likely results for planetary systems not modeled in this 
section. 

An important question for the observing teams is how long they must continue frequent 
monitoring in order to have a good chance of detecting a repetition. Figure 9 indicates that 



^°These rates are computed using the normalization described in §3.2 and §4.2. That is, they are 100 
times the number of events of the type described, divided by the number of events in which the central star 
served as the lens in an isolated event, with the peak magnification of the stellar event achieving a value 
greater than 1.34. 
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Double— Peak Events: = Gaussian, a = uniform, detect 
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Fig. 8. — Distributions of encounter durations and of and wait times between encounters for 
events with one repetition, using the detection criteria of set A. Durations: For the solar 
system (top left), Saturn dominates the left- most peak, with Uranus and Neptune appearing 
as a shoulder to the left, at approximately 0.7 days. Overlap doubles are shown with a 
dotted line: overlap encounters involving both Jupiter and the Sun are represented in the 
right-most peak; the other encounter that occurred as part of these repeating events was 
with either Saturn, Uranus, or Neptune, and produced the peak on the left. Wait times: 
For the solar system, the predominant peak is produced by doubles involving the Sun and 
Saturn; Uranus and Neptune produce a broad shoulder to the right, at around 10^ days. 
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Double-Peak Events: Vj = Gaussian, a = uniform, detect = A 
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Fig. 9. — Each panel shows the integrated area under the distribution to which it corresponds 
in the previous figure. Encounter durations are shown on the left; wait times between 
encounters are shown on the right. Roughly 75% of all encounters, and half of the gaps 
between encounters (75% for the solar system), last less than 100 days. 
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Double-Peak Events: Vt = Gaussian, a = uniform, detect = C 
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Fig. 10. — Distributions of durations of and wait times between events with one repetition, 
using the detection criteria of set C. Most of the structure visible in Figure 8, where the 
detection criteria of set A were used, can also be seen here. The overall event rates have 
increased, and the gaps in the centers of the distributions have filled in. A new peak 
corresponding to Pluto has appeared in the duration distribution for the solar system (top 
left); these events are so short that, although we may occasionally catch one magnified point 
on their light curves, we are unlikely to be able to follow them. 
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1/2 — 3/4 of all events with one repetition will begin the repetition within 100 days of the 
time the magnification associated with the first encounter dips below 1.06. In addition, a 
large fraction (~ 75%) of overlap events may be identifiable. In our simulations, overlap 
events form a significant subset (1 — 4%) of all events for all of the models (solar system 
through to power-of-2) and for all detection criteria (sets A through C). This may be an 
artifact however, because in each case the model system happened to have a planet located 
just outside the zone for resonant lensing. Note though, that if overlap events are not seen 
at roughly this level, there are hkely to be few resonant events, since the planets that serve 
as lenses in overlap events are the ones that can be brought into the zone for resonant 
lensing as the orientation of the planetary system changes. 

5. Microlensing and the Search for Extraterrestrial Life 

Recent and ongoing advances in technology have led to the discovery of 
extrasolar planets (see references hsted in the Encyclopedia of Extrasolar Planets, 
www.obspm.fr/darc/planets/encycl.html), and promise the discovery and even the imaging 
of additional planets (Angel & Woolf 1997; Praclas & Shelton 1997; Labeyrie 1996; Brown 
1996). These developments excite the imagination because they seem to bring us closer 
to the possible discovery of extraterrestrial life. It is therefore interesting to ask whether 
microlensing searches are likely to find planets on which life could thrive. 

We do not yet have a clear enough understanding of the nature of life to definitively 
answer this question, because the range of physical conditions compatible with life may 
well be wider than our limited experience would at first suggest. It has been proposed, for 
example, that life may exist on the outer planets of our own solar system and/or on their 
moons. (See, e.g., Reynolds et al. 1983; Raulin et al. 1992; Sagan, Thompson, & Khare 
1992; WiUiams, Kasting, & Wade 1997, McCord et al. 1997.) It has even been postulated 
that life may exist in non-planetary environments, including the interiors of stars and 
molecular clouds (see, e.g., Feinberg & Shapiro 1980). Thus, any planet, however close to 
or far from its star, and whatever the nature of the star, can be considered as a possible 
harbor for life. 

Nevertheless, in the absence of real information on the existence of life away from our 
own planet, one question is clearly interesting: will microlensing find evidence of planets 
similar to Earth? We must of course define what we mean by "similar to Earth" . If we 
mean that there is a chance that chemical processes necessary for Earth-like life could 
occur, then we want to consider planets which can have similar surface and atmospheric 
make-up, and similar amounts of energy available to fuel the necessary chemical processes. 
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Thus we may ask: by how much can we change the mass of the Earth, Mq, and its distance 
from the sun, a® (i.e., the quantities microlensing events can measure) and still maintain 
conditions that seem able to support life forms with chemistry similar to that of Earth life? 

The range of planetary masses and distances from a solar-type star compatible with 
Earth-like conditions, particularly the presence of liquid water, has been dubbed the 
"Goldilocks Problem", and has been studied by many researchers. (See, e.g., Rampino & 
Caldiera 1994, for an overview.) Here we make some general observations. 

The primary requirement on a is that the incident flux of radiation from the star should 
be comparable to the flux received by the Earth from the Sun. That is, T jT^^ should not 
be too different from unity. 

Uo; \M)\Ds) \x{x-\)) ^ ^ 

Here we have assumed that there is a planet located a distance [iRe from the central star; 
the mass and luminosity of the central star are M and L, respectively. As usual, Dg (Dl) 
is the distance to the lensed source (lens), and x — Dl/Ds- Given that (1) the conditions 
that lead to life may be flexible, (2) the effects of radiation incident from the star are likely 
to be strongly influenced (either enhanced or diminished) by the planet's atmosphere, and 
(3) internal heating from geological processes or radioactive materials may be important, it 
is not clear how large a range of values of J-'/J-'^ may be compatible with the development 
of life. We will therefore simply use T jTc^ = 1 as a guideline. 

Figure 11 shows the relationship between M and x for those systems that satisfy the 
relationship T jT;^ = 1. We have assumed that L = M^'^ in the upper plot, and L = 10 M^'^ 
in the bottom plot; the former would be appropriate for main-sequence stars, while the 
latter would be appropriate for slightly evolved stars. As Eq. 18 and Figure 11 make clear, 
if the planets we will discover via microlensing are to have an incident flux comparable to 
the flux incident on Earth, their stars will generally (although not necessarily) be more 
luminous than our sun. This has two obvious implications. The first is that the length of 
time during which the planet would have this flux incident will be shorter than the time to 
date that the Earth has had roughly this flux incident. This is because the system's star 
may need to be more massive than the Sun, or even slightly evolved. The time elapsed from 
the formation of the star until the present time could range from less than 0.1 the present 
age of the sun to times comparable to the sun's main-sequence lifetime. We do not know 
how long it takes for complex life forms to develop, but it may be that the process is fast 
enough that intelligent life can develop and thrive during a time significantly shorter than 
the present age of the Sun. Indeed, it is likely that a long sequence of independent processes 
must occur in order for intelligent life to develop; thus, the probability distribution may 
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Fig. 11. — Plotted is M vs x, where M is the mass of the central star of the lens system 
and X — Dl/Ds- We have assumed that there is a planet located at a — /iRe, with /i — 1.5, 
and have imposed the condition: T jT® — 1. Upper curve in each panel: L — M^-^. Lower 

curve in each panel: L = lOM^'^. For the Bulge [Magellanic Clouds, M31], we have assumed 
Ds — 10 kpc [60 kpc, 700 kpc]. The plot in the upper left panel is for the Galactic Bulge; the 
upper right panel is a blow-up of the region in x that corresponds to the lens, as well as the 
lensed source, being located in the Bulge. Horizontal dotted lines are drawn at M = Mq and 
M = 2 Mq The left-right pairs in the middle and bottom panels show the same quantities 
for the Magellanic Clouds and M31, respectively. 
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well be log normal, and the likelihood of intelligent life developing in times much shorter 
than the time apparently taken on Earth may be significant. 

The second implication is that, since the central star must be fairly luminous, it may 
contribute a non-neghgible fraction of the hght incident along the hne of sight to the lensed 
source; i.e., the light we receive may be strongly blended. We return to this issue in the 
next section, but note here that it is potentially important, because when the planetary 
system's star emits enough light to influence the lensing light curve in a detectable way, we 
may be able to learn more about the planetary system. 

We may also argue, as follows, that planets hkely to be deemed "Earth-hke" have 
masses within a factor of ~ 15 of the mass of the Earth. Assuming that we would like a 
rocky surface, we also assume that the planet's average density should be similar to that of 
Earth. This means that the acceleration due to gravity, scales as the cubed root of the 
planet's mass, m. Thus, g will be within a factor of 2.5 of g^) (= 9.8 m/s^) if m is within a 
factor of 15 of m®. Increasing or decreasing the value of g will lead to different atmospheric 
contents; for any given atmospheric temperature, there is a lower limit to m (hence g), 
below which an atmosphere will not be retained. This means that, excluding finite-source 
size effects (which can significantly increase the duration of some events), event durations 
for Earth-like planets are likely to range from ~ 2 to ~ 32 hours. 

We should also consider the possibility that microlensing (whether by resonant or wide 
planets) is most likely to discover the outer planets in a system that may contain closer 
planets experiencing Earth-like conditions. In this case, the microlensing events serve as 
beacons directing us to the planetary system the outer planets inhabit. Whether the planets 
that serve as lenses themselves harbor life, or whether they merely serve as bellwethers, it 
is important to consider whether additional observations can help us to learn more about 
the planetary system. 

6. Blending 

If the central star of a planetary system that serves as a lens is fairly luminous, then 
its light will blend with that from the lensed star. Let / represent the fraction of the 
baseline flux contributed by the lensed star. We expect that in most cases, / is a function of 
wavelength. If studies of the lensing light curve, or spectra taken during the event provide 
evidence that light from the lensed star is blended with light from other sources, we can 
hope to measure the value of / in each of several wavebands. (See, e.g., DiStefano & Esin 
1995.) This information, combined with our study of the baseline flux, allows us to obtain a 
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spectral type for the lensed source and also to study the spectral character of the unlensed 
light. If we are then able to determine that the unlensed light, or at least a quantifiable 
portion of it, emanates from the star in the lens system, we may then be able to determine 
the spectral type and possibly the mass of that star as well. If we are able to repeat this 
process for a set of planetary-system events, we will thus be able to gather statistics about 
the types of stars which, in the distant reaches of our Galaxy, or in other galaxies, have 
planetary systems. When the hght curve allows us to determine the mass ratio q between 
at least one planet and the central star, then the mass of that planet is also determined. 
Note that the process described above can be carried out even if the track of the source 
passes close to a planet, but not close to the central star. 

When the central star of the planetary system is luminous, particularly if its flux comes 
close to satisfying the criterion studied in the last section {T jT^ ~ 1), then / can indeed 
be small enough for the effects of blending to be measurable. If, for example, the apparent 
V magnitude of the combined light coming along the line of sight from a lensing event is 
[19.5, 17.0, 15.7], then / < 0.1 if the lens is located in the Bulge and is a main-sequence star 
of mass roughly equal to [1, 2, 3] Mq. A small value of / can allow us to reliably determine 
the effects of blending and to thereby learn more about the lens. Detection of the event 
can, however, become somewhat more difficult. 



When light from the lensed source is blended with light from other sources, the 
observed magnification, Aq^^, is smaller than the true magnification, A. 



Thus, in order for a light curve perturbation to be brought above the detection limit, A 
must be larger than it would otherwise have to be, the projected distance between the 
source and lens must be smaller, and the event will consequently appear to have a shorter 
duration. To describe this effect systematically, DiStefano & Esin (1995) introduced the 
"blended Einstein radius", RE,h- As before, let be the minimum peak magnification 
needed for event detection. Here, however, we consistently define the duration of each event 
to be time during which the magnification was greater than The expression for the 

blended Einstein radius is then 



6.1. The Effects of Blending on Event Detection 



(A„6,-1) = /(A-1) 



(18) 



-Re,6 — Re^/^ 



{A„iin 1) + / 



- 1 



(19) 
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Figure 12 illustrates the influence of blending on event duration as a function of both 
Ajnin and /. For a given value of /, the event duration is longer if is smaller. Thus, 
increasing the photometric sensitivity to accomodate smaller values of should increase 
the detection rate. Note that if Amin is 1.06, then even if / = 0.18, the event duration 
(which would have been 2 te), is reduced by only a factor of 2 to 1 te- Thus, while blending 
does tend to decrease the event duration, making more frequent monitoring desirable, we 
can expect to be able to detect a large majority of events by using sensitive photometry. We 
note that, if Amin could be reduced to 1.02-a formidable task for a large-scale monitoring 
program-then for / = 0.03, an event that would have lasted for 3te will have an observed 
duration of 1 r^;. |^ 



6.2. Studying Blending 

As long as the presence of a blend of light from sources other than the lensed star does 
not prevent us from detecting or identifying an event, blending can be a great boon to the 
analysis of microlensing events. Using blending to full effect requires several steps. 

(1) Spectra: Spectra taken during the event can help us to determine the spectral type of 
the lensed source. During the event, the lensed source contributes a fraction, fevit), of the 
incident light. 

Thus, spectra taken at peak can be compared with spectra taken at baseline (and possibly 
at other times during the event) to determine the spectral type of the lensed star, and to 
thereby better quantify the contributions of any unlensed source of light as a function of 
wavelength. 

(2) Light Curve Fitting: If the light curve is observed in several wavebands, then the flt to 
the light curve will determine / in each waveband separately. This information, combined 
with the baseline magnitudes, also provides a way of determining the spectral type of the 
lensed star, and therefore the spectral characteristics of the unlensed light. Information 
about blending drawn from studying the light curve, should confirm and complement 
information drawn from spectral studies. 



^^Wambsganss (1997) has illustrated that, when there is a planet in the zone for resonant lensing, evidence 
of the planet's presence can be exhibited relatively early during the rise of the underlying stellar- lens event. 
In such a case, a smaller value of Amin will also help us to discover resonant events that happen to be subject 
to blending. 
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Attenuation of Event Duration 
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Fig. 12. — Log[f] vs Log[Aj„j„]. Source tracks going straight through an Einstein diameter 

have been considered, te is the time needed for the source to cross a distance equal to 
the Einstein diameter. The observed event duration, r, is defined to be the time during 
which A > Afnin- Because each event was subject to blending (/ < 1), r is attenuated 
relative to the event duration without blending. Note that, for A„jj„ < 1.34, the duration 
of the observed event would be greater than te, were there no blending. Each curve shown 
corresponds to a fixed ratio t/te- The horizontal fines correspond to fv — 0.5,0.1,0.01, 
proceeding from top to bottom. For reference, we note that if the lens is a main-sequence 
star with mass [1, 2, 3] M©, then fv — 0.5 if the total apparent V magnitude along the fine 
of sight is [20.2,17.6,16.4]; similarly, fv = 0.01, for My = [19.4,16.9,15.6]. Vertical lines 
correspond to Amin = 1.02, 1.04, 1.06; these provide a comparison which demonstrates the 
improvement in event detection relative to the case A — 1.34. 
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(3) High-Resolution FoUow-Up: The determination that there is hght incident from an 
unlensed source does not necessarily indicate that the extra hght is emanating from the 
lens itself. In fact, the fields studied by the monitoring teams are crowded, so there is a 
high probability that another star, independent of both the lens and lensed source, may 
be contributing to the incident flux. High-resolution observations can determine whether 
the "extra" light is emanating from one or more stars, distinct from the lensed source, that 
were not resolved by the monitoring teams. The density in the Bulge, for example, is small 
enough that it is unlikely that independent stars located along the direction to the event 
(apart from the lens), will not be resolved by HST. Because the angular separation between 
the lens and source is on the order of milliarcseconds, we cannot, however, resolve the lens 
and source immediately after the event, even if the lens system contains a luminous star. 
Thus, the light from each star resolved by high-resolution follow-up can be subtracted (as 
a function of wavelength) from the baseline flux to determine how much light, in each of 
several wavebands, may be emanating from the lens system itself. After a time determined 
by Vt, we should be able to resolve the lens system and test the hypothesis that the lens 
provides the light unaccounted for by other resolved stars; the wait-time is typically on 
the order of decades. Until the lens and source are resolved, we cannot be certain that the 
additional hght does come from the lens system. It is after all possible that yet another 
star has a small angular separation from both lens and lensed source; the likelihood of such 
a coincidence can be evaluated. In many individual cases, and certainly for a statistical 
sample of events, we expect that the most likely alternative is that the additional light 
emanates from the lens system. Spectral and high-resolution measurements may provide 
enough information to determine the spectral type of the central star of the planetary 
system. Although the distance from the lens will not be known, we will, particularly with 
the input from the lensing event itself, be able to place limits on its mass. (See Udalski et 
al. [1994] for a discussion of this problem for the OGLE 7 event.) The mass of the central 
star, together with the determination of q for any planet lens, allows us to compute the 
planet's mass. 

(4) High-Spectral-Resolution FoUow-Up: We note that Keck HIRES observations have 
produced a high-resolution spectrum for a star (with My — 17) in Baade's window, allowing 
a detailed abundance analysis to be carried out (Castro et al. 1996). This would also be 
possible for stars in the planetary systems we discover via microlensing. Ideally, one would 
like to go further, and to achieve spectral resolution good enough to detect evidence of 
low-velocity motion, which may indicate the presence of planets closer to the central star 
than the planets discovered via microlensing. Extrapolating from the HIRES work and 
from more general considerations about radial velocity measurements (Latham 1996), we 
estimate that the photon flux we are likely to receive from the lens system is too small 
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(generally by more than 2 orders of magnitude) for the present generation of telescopes 
and detectors to be effective. Thus, although possible in principle, such fine radial velocity 
measurements are presently beyond our reach for lens-system stars in the Bulge. 



7. Finite-Source-Size Effects 

Finite-source-size effects are known to be important when attempting to detect 
evidence of a low-mass (e.g.. Earth- mass) planet located in the zone for resonant lensing; 
in this case, it is the comparison between the size of the caustic structures and the radius 
of the source as projected onto the lens plane, that is important. The caustic structure, 
however, shrinks and therefore plays a less important role as the separation between the 
star and planet increases. Thus, for planets in wide orbits, the important comparison is 
between the size of the planet's Einstein radius and the size of the source as projected onto 
the lens plane. Until now we have assumed that the size of the lensed source can be 
neglected in comparison to the size of the Einstein ring. For lensing by planetary systems, 
this is not always a good assumption. 

Let r represent the ratio between the radius of the source star and the Einstein radius 
of the lens, r = Rs/RE,i, with Rs the radius of the lensed source as projected onto the lens 
plane. 
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where M is the mass of the lens. If the lens has the mass of Jupiter and if the source and 
lens are both in the Bulge {Ds = 10 kpc, x = 0.9), then r ~ 0.05. Witt & Mao (1992) 
found that the finite size of the source begins to influence the shape of the light curve 
when r ~ 0.1. Since typical source stars are likely to be larger than Rq, finite source size 
effects may play some role even when a planet as massive as Jupiter serves as a lens. The 
significance of their role increases as the mass of the planet lens decreases. 



7.1. Finite Source Size and Detectability 

The light curve changes in two ways when the size of the source cannot be neglected. 
First, whereas the magnification of a point source lensed by a point mass can become 



^^The caustic structure may play a more important role in a small fraction of overlap events than it does 
for wide-orbit events generally. 
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arbitrarily large when the distance between the lens and projected source position becomes 
arbitrarily small, the magnification has an upper limit when the source size is finite. Thus, 
the peak of the light curve becomes attenuated. The larger the source size relative to 
the Einstein radius, the more pronounced the attenuation. The maximum magnification 
possible for a disk of constant surface brightness, for example, is 



This formula indicates that the size of the source must be fairly large in order for the peak 
magnification to be brought below the level of detectability, particularly if the photometric 
sensitivity is good. For example, a magnification of 1.34 can be achieved even if the size 
of the source is ~ 2.25i?E, and a magnification of 1.02 can be achieved even if the size of 
the source is ~ IORe- Note that a projected source size of 10-R^; corresponds to a ~ 200Rq 
{20Rq) stellar radius when the lens is of Jupiter-mass (Earth-mass). (We have assumed 
X > 0.9.) Thus, in spite of the decrease in peak magnification, a goodly fraction of stars in 
the source population can, when are lensed by planets in wide orbits, produce events that 
should be detectable, at least with good photometry. 

The second effect of finite source size on the light curve is to broaden the width of the 
perturbed region of the light curve, since some portions of the source may be significantly 
magnified well before the center of the source achieves its closest approach to the lens. The 
light curves can become almost fiat-topped. Thus, magnifications close to the peak value 
may be sustained during the time it takes the source to travel several Einstein diameters, 
and the event appears to last longer than it would have had the lensed source been a point 
source. This can be a boon to the detectability of events due to planet-mass lenses, since 
one of the greatest barriers to detection is the short lifetime of the event. 

Figure 13 displays a sequence of light curves for the case Rs — 5Re- When comparing 
the different light curves, note that the peak magnification does not change much, even 
as the distance of closest approach, Dj^in: changes from to 4: Re- If the photometric 
sensitivity of the observations was such that a events with peak magnification above 1.06 
could be detected, then all of these events would be detected. The primary difference as 
Dmin increases, is that the time duration of the observed event decreases. Even so, for Dmin 
as large as 4: Re, the observed duration is still roughly equal to the time a point-source 
light curve would remain above A = 1.06. The relationship between -Dmin and Rs is 
explored more systematically in Figure 13. The implications are that: (1) depending on the 
distribution of source sizes, the detection probabilities can be increased by finite-source-size 
effects, and (2) events that can be detected can also generally be monitored for longer 
periods of time. 
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Fig. 13. — In each panel, the projected size of the source in the lens plane is 5RE,i- The 
distance of closest approach, Dmin, decreases from ARs^i in the upper panel to the in the 
lower panel. Note that the peak magnification is larger than 1.06 in all cases. Furthermore, 
if the measured event duration is the time during which the magnification is larger than 
1.06, then these events last up to 2.5 (bottom panel) times as long as they would have, had 
the lensed source been a point source. 
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Fig. 14. — The source radius, Rs, is plotted against -Dmm, the distance of closest approach. 
Both quantities are expressed in units of the Einstein ring radius. Each point corresponds to 
an event in which the peak magnification was roughly equal to the value listed on the upper 
right of each curve. For each value of Rs, only the points with the largest values of 
leading to the listed magnification are shown. We have assumed that the disk of the lensed 
source has constant surface brightness. The value of Dmin for a point source is shown, for 



1.02, 1.06, and 2.0 [vertical lines]. 
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Fig. 15. — Results of a Monte Carlo simulation that includes finite-source-size effects for a 
lensing system with 9 Earth-mass planets whose separations from the solar-mass central star 
are the same as those of the planets in our solar system. We have assumed that the source 
radii are distributed uniformly from to 10 Re- In addition, wc assume that lensing by 
Earth-mass planets cannot be detected if Rs > 7 Re- Thus, roughly 43% of the sources 
could not produce identifiable lensing events. For the remainder of the sources, we used 
the approximation that D^i„ ~ Rs- (See Figure 14.) In these simulations we used: Vf — 
Gaussian, a — uniform, detect = C. For comparison, the distributions for the equivalent 
simulation of lensing by the model of the solar system used in §4, with the sources taken to 
be point-like, are shown by the dotted curves. 
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Gould (1994) was able to find point-source fits to light curves associated with 
large sources. This degeneracy of the light curve shape, however, need not prevent the 
identification of individual events as having been influenced by finite-source-size effects, 
since spectral information can break the degeneracy. Real stellar disks exhibit brightness 
profiles and spectra that have spatial structure. When a star for which r is not negligible is 
lensed, the observed spectrum will therefore be time-dependent, and can be used to study 
the lensed star. (See, e.g., Witt & Mao 1994, Witt 1995, Gould 1994, Loeb & Sasselov 
1995, Simmons, WiUis & Newsam 1995, Sasselov 1997, Heyrovsky & Loeb 1997, and, 
for an observation, Alcock et al. 1997.) This time dependence also helps to confirm the 
interpretation of the event as due to microlensing of a source whose size has influenced 
the light curve. Indeed, when the source is large relative to the Einsein ring, it is as if 
a magnifying glass were scanning the face of the lensed star, thus providing a good deal 
of information about the source star and about the event. Spectral studies, even during 
short-duration events can therefore be valuable in breaking the degeneracy for individual 
events. Even if spectral information is not available, however, statistical analysis of the 
ensemble of events would be able to indicate that flnite-source-size effects had played 
a role. In particular, the relationship between peak magnitude and duration would be 
different than it would be for a set of true Paczyiiski fight curves; for large sources, the peak 
magnification is closely related to source size, while the duration continues to have a more 
direct relationship to the distance of closest approach. 

To better quantify the influence of flnite source size, we have carried out simulations 
using the solar system model and a power-of-3 model to study the detectability of 
Earth-mass planets. These results are summarized in Table 4 (§10) and are illustrated for 
the solar system model in Figure 15. (The input assumptions are described in the caption 
to Figure 15.) Because their Einstein rings are smaller, Earth-mass planets are expected 
to yield a smaller event rate. The rate would be smaller than that due to Jupiter-mass 
planets by a factor of roughly 18. Yet, in the power-of-3 model, where all of the planets in 
our original model (§4) had a mass equal to that of Jupiter, the computed attenuation of 
the Earth- mass system relative to the Jupiter- mass system is only a factor of ~ 3. Even 
though flnite-source-size effects made event detection impossible when the Earth-mass 
planets lensed roughly 43% of the sources, the fact that the detection rate was higher 
for the remainder of the sources provided a net increase in the detection rate over what 
might otherwise have been expected. The relative increase is even more pronounced for the 
solar-system model, since Saturn, Uranus, Neptune, and Pluto are responsible for many 
of the isolated events when a model of the solar system (§4) serves as a lens, and these 
each have a smaller Einstein radius than Jupiter. Although more realistic stellar surface 
brightness profiles should be used, the basic result that Earth-mass planets are detectable 
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should be robust. Also robust is the increased event duration for a range of values of Rs- 
For a fixed value of Amini this corresponds to a larger value of the effective width, Wi. 
and therefore a larger value of n = Wi/RE.i- This can increase the detection probabilities 
above what is expected if the source is point-like, as predicted and as demonstrated by the 
simulations. 



When spectroscopic studies of the lensed source can establish its real physical size, 
finite-source-size effects may then allow us to derive the true size of the Einstein ring. 
This means that, if the value of r can be measured from the data, Vt can be determined. 
Furthermore the degeneracy in the lens mass associated with the Paczyhski light curve is 
partially broken, since knowing the value of Re, and, through spectral studies, the value of 
Ds, allows us to infer 



where C is a constant whose value is measured. Given the distribution of likely values of x 
(from 0.9 to 0.99, for example, in the Bulge) this generally constrains M to within a factor 
of ~ 5. Thus, the hypothesis that any given event is due to lensing by a planetary-mass can 
be meaningfully tested. 



Expanding the microlensing searches for planets to those in wide orbits allows the 
presence of each wide-orbit planet to contribute to the probabihty of detection, and also 
allows for possibly complex signals, containing information about the masses and spatial 
distribution of several planets in the planetary system. It is therefore interesting to consider 
that the complexity of planetary systems extends beyond the presence of multiple planets, 
to the possible existence of moons revolving about the planets, and also of belts or rings 
of space debris, such as asteroids and comets. What is the probability of detecting these 
features through microlensing studies? 



7.2. Planeteiry Mcisses and Finite-Source-Size Effects 
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8. Moons, Asteroids, and Comets 



8.1. Moons 



The microlensing detection of a moon about a planet requires that (1) the planet about 
which it revolves either be in the zone for resonant lensing or in a wide orbit, (2) the moon 
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be in either a resonant or wide orbit about the planet, (3) the duration of the perturbation 
due to the moon is long-hved enough to be detectable, and (4) finite-source-size effects do 
not wash out the signature of the moon's presence. 

8.1.1. Planets in the Zone for Resonant Lensing 

If the planet is in the zone for resonant lensing, then the presence of the moon will 
typically complicate the caustic structure and the associated pattern of the light curve. 
For planets in wide orbits, the presence of a moon in the zone for resonant lensing about 
the planet could mean that the light curve associated with lensing by the planet is more 
likely to exhibit the behavior associated with caustic crossings. The size of the caustic 
structures due to the presence of a moon will generally be small; consequently the associated 
perturbations may be difficult to detect. 

8.1.2. Wide Orbits 

If the moon is in a wide orbit about a planet which is itself in a wide orbit about its 
star, it is relatively straightforward to determine what the signature of the moon's presence 
is, and to assess its detectability. Going through the 4 criteria listed above, we note that 
the first two conditions are demonstrated to be feasible by the example of our Solar System. 
Were an identical system to be placed in the Bulge, Saturn and Neptune would both be in 
wide orbits, with projected separations significantly larger than the sun's Einstein radius 
most of the time and for most system orientations. Pheobe orbits Saturn in an orbit that 
is wide with respect to Saturn's Einstein radius; similarly, Nereid orbits Neptune in a wide 
orbit (Norton 1989). With regard to criterion 3, if we would like an event due to lensing by 
a moon to last at least 1 hour, then, assuming that an event due to the central star lasts 
100 days, the mass ratio between the moon and the star can be as small as 1.7 x 10~^. If the 
central star has a mass equal to that of the sun, the smallest detectable moon could have a 
mass as small as rrimin = 3.5 x 10^^ gm. This is roughly 5 times the mass of our own moon. 

Criterion 4, the limit placed by finite-source-size, is also important. If we impose the 
condition that r have a maximum value of [10, 5, 1], then we find that the minimum moon 
mass has a value of [5 x 10^^ gm, 2 x 10^^ gm, 5 x 10^'' gm]. These simple calculations make 
it clear that finite-source-size effects will be important in the microlensing study of moons. 
Those moons that do give rise to detectable microlensing events may yield light curves 
whose shape is significantly infiuenced by finite-source-size effects. In these cases, the effects 
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can be beneficial, since the event durations may be prolonged by finite-source-size effects, 
lasting several hours. 

We note that, if its density were 6 gm/cm^, a moon of mass equal to 3.5 x 10^^ gm 
would have a radius equal to ~ 2400 km. In fact, the density of the wide moons in our 
own solar system is estimated to be smaller than 6 gm/cm^, while their physical radii are 
smaller than 1000 km (Norton 1989). It is therefore unlikely that microlensing observations 

of our solar system would reveal the presence of its many moons. There does not seem 
to be a fundamental reason, however, why other systems might not have moons that can 
be detected via microlensing. Thus, it may be productive to institute even more frequent 
monitoring just after a short-duration event; if the event is due to a planet, careful frequent 
monitoring after it has ceased would improve our ability to detect any associated moons. 



8.2. Asteroids and Comets 

It is also interesting to consider the detectability of a belt of debris, such as the solar 
system's asteroid belt, or a cloud of comets. Each isolated member of such a belt is likely 
to be much too small to be detected through its action as an individual lens. It may be 
possible, however, that the projected mass surface density is close to the critical value, 
acrit- In this case, the track of the source could pass across a caustic structure, giving rise 
to variability in the light curve. (See also Keeton & Kochanek 1997.) The time scale of the 
perturbations would be set by the time required for the track of the source to pass through 
the belt. 

The critical density is 

= 4Gj/J(l-x) = " B) (^) (tJT^) 

This surface mass density would not be achieved by the belts of asteroids and comets 
associated with our own solar system, unless Ds is much larger than the distance to the 
Bulge, and the value of x is neither too large nor too small. 

Consider a belt of debris of total mass 10'^" gm, with a projected extent of 0.01 AU 
by 1 AU. The surface density would be ~ 4.4 x 10^ gm/cm^. With Ds equal to 10 (60) 
kpc and x = 0.5, we derive (Tcrit =~ 4.4 x 10^ (~ 7.3 x 10^ gm/cm^). Thus, it is not 
inconceivable that, early in the evolution of a planetary system, acrit can be achieved, 
leading to interesting signals, particularly for favorable system orientation relative to the 
line of sight. 
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9. Detection Strategies 

The strategies best suited to the discovery of wide planetary systems have much in 
common with those best suited to the discovery of planets in resonant orbits. In both cases 
it is important to have frequent time samphng of events or event features whose duration 
may be on the order of hours. In both cases, good photometry is important if we are to 
carry out meaningful tests of the planetary-system lens hypothesis. Thus, since a detailed 
search mechanism is already in place to discover planets in the zone for resonant lensing, 
it is useful to ask how effective that mechanism will be at discovering wide planets, and 
whether slight modifications of it will make it more effective. 

9.1. Past and Ongoing Observations 

9.1.1. Short-duration events 

The monitoring teams have discovered some short-duration events, with the shortest 
event on record lasting approximately 2 days. To better understand what is known so far, we 
have considered events along the direction to the Bulge, since the MACHO team has already 
accumulated a store of over 150 Bulge events about which some information is publicly 
available. Specifically, we have used their paper on 45 Bulge events monitored in 1993 
(Alcock et al. 1997a), and their "alert" web pages ( [ittp: / / darkstar .astro. Washington. edu|) , 
which list approximate durations for many of the 136 Bulge events observed during 1995, 
1996, and 1997. Altogether we found 148 events (a) which are apparently due to lensing by 
a point-mass, and (b) for which approximate durations are available. (We note however, 
that the durations posted on the web site come from fits that may be refined in the future.) 
The distribution of event durations is plotted in Figure 16. Note the appearance of the 
short- duration peak. In addition, at least one short-duration event is not included in this 
graph; 96-BLG-9 is simply listed as a "short-time-scale event" , with no duration given. 

In light of Figure 16 we may ask whether the MACHO team has already begun to 
discover concrete evidence for a signature due to planets in the Bulge. This may be the 
case. (In fact Bennett et al. [1996] have conjectured that one short event may indeed be 
due to a planet in a wide orbit .) We emphasize, however, that a good deal of further work 
would be required to clearly establish it. In fact, the remainder of this section and the next 
are devoted to detailed discussions of productive lines for further research. One certain 
conclusion we can draw from the MACHO Bulge data is that the team has proved that it 
has the ability to detect short- duration events. This is encouraging, since the detection 
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Fig. 16. — Distribution of event durations for the MACHO team's set of Bulge events. 
Included is data from 1993 (Alcock et al 1997a), and the "alert" web page data for events 
which occurred in 1995, 1996, and 1997. All 148 events which were not hsted as binary 
events, and for which estimates of event durations had been, made are included. 
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strategies presently used are not optimized for the discovery of short time-scale events. 
Thus, the discoveries to date indicate that the number of short duration events is large 
enough that a more concentrated effort to find them is likely to bear fruit. 

In addition to being able to discover short-duration events, the MACHO and EROS 
teams have also been able to place upper limits on the numbers of such events along the 
direction of the LMC (Renault et al. 1997, Ansari it et al. 1996; Alcock et al. 1996). 

The follow-up teams can only study short-duration events discovered by the monitoring 
teams if the events last more than a day. In fact, the time presently taken to call an alert 
may vary from 1 day up to ~ 5 days (Axelrod 1997). Thus, the chances that the follow-up 
teams will be able to catch lensing events due to planets with mass less than that of Jupiter 
are small. They may, however, be able to follow-up on events due to Jupiter-mass planets 
and larger, or on events due to less massive planets if the durations of the events are 
extended by finite-source-size effects. 



9.1.2. Repeating events 

The early folklore on repeating events said that they should be discarded: the 
micro lensing of light from a specific star is such a rare event that, if a microlensing-like 
signal is seen twice when monitoring the same star, the repetition is almost surely a sign 
that the observed variation is due to something other than microlensing. Although it has 
been pointed out that binary sources (Griest & Hu 1992) and binary lenses (Di Stefano 
& Mao) can each lead to microlensing events that repeat, it is nevertheless the case that 
one of the cuts used to eliminate lensing candidates from further consideration is based 
on whether the event appears to repeat. As a result, the monitoring teams may not 



This is an appropriate place to note that, while our choice of 1 day for the minimum time-duration 
necessary to detect a first encounter (detection criteria A and B) is perhaps too optimistic to describe the 
present search strategy, the fact that the monitoring teams have been able to identify events that last for 
~ 2 days indicates that it may not be overly optimistic. It can clearly be achieved in some fields even 
at present; were more emphasis to be placed on the discovery of short-time-scale events, then it might be 
achievable for a larger portion of the data sets. In addition, as indicated earlier, even if frequent monitoring 
of a short-duration event does not begin in time to provide much more information about its light curve, it 
could nevertheless play a crucial role in discovering any later repetitions. 

^"^Of course the complication of ensuring that a repetition is not due to stellar variability is serious, since 
variability is much more common than microlensing. Nevertheless, if a star that has never been observed to 
vary from baseline produces two deviations well-fit by microlensing models, it is certainly worth exploring 
the possibility that the deviations are indeed due to microlensing. In some cases blending or finite-source-size 
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pursue the implications of evidence of repeating events in their data sets. We note that one 
"repeater" is mentioned on the MACHO alert web site, but the posted information does 
not indicate whether microlensing seems a likely explanation for either perturbation. 

The follow-up teams do not, at present, continue to follow events once the detected 
flux has returned to its baseline value. 



9.2. Useful Modifications 

9.2.1. Improved Measurements for both Resonant and Wide Planets 

Since detection of blending and finite-source-size effects could possibly clinch the 
planet-microlensing interpretation, it is worthwhile to search for evidence of these effects 
in every event. The analysis of light curve shape can play an important role; this analysis 
can be done after the event has ceased. During the event, it would be valuable to obtain 
spectra, particularly near the time of peak magnification. 

9.2.2. Improved Measurements for Wide Planets 

The main difference between the resonant and wide cases is whether or not there is 
a trigger that indicates that frequent monitoring should begin. In the resonant case, an 
ongoing stellar lens event provides a convenient trigger. As soon as any microlensing event 
is identified, intensive monitoring that lasts for the duration of the event can begin. In the 
wide case, many of the events that could be better understood through intensive monitoring 
have no convenient trigger. This is particularly true for any event in which a planet is the 
first or only lens. Thus, in order to optimize our ability to detect planet-lens events, regular 
frequent monitoring is called for. The optimal frequency however, depends on the 
minimum peak value of the magnification for which an event can be reliably identified. This 
dependence is best expressed in terms of n = Wi/RE^i- If k is the number of repetitions 
exhibited by the light curve, the detection rate scales as n^+^. Not only arc more events 
discovered when n is greater than unity, but the duration of detected events is also longer, 
making less frequent monitoring possible. Note that finite-source-size effects can increase 
the value of Wi (hence n), even for a fixed value of Amin- 



effects may make it possible to carry out further tests of the conjecture that multiple disturbances are due 
to microlensing. 
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The bottom line is therefore that good photometry is the key to increasing the 
detection rate of planets in wide orbits. Monitoring that is consistently more frequent 
than that presently carried out by the monitoring teams is also important, although hourly 
monitoring may not be needed even for the detection of Earth-mass planets. 

The design of the observing programs of the future should take these requirements 
seriously. There are likely to be several different ways to accomplish them. For example, 
a system of monitoring could achieve a detection rate ~ 10 times as high as the present 
system does by (1) using photometry as good as or better than the present follow-up teams, 
(2) using pixel techniques to identify lensing of stars that are otherwise (i.e., at baseline) 
below the detection limit, and (3) monitoring frequently. This would mean discovering 
C(IOOO) Bulge events per year. Not only would the numbers of events be high, allowing us 
to detect some very low- probability events (such as multiple repetitions), but the quality 
of the data would be so high that we would understand each event and the population 
of lenses and sources much better than we do at present. Another approach might be to 
have a world-wide network of telescopes, each taking deep images of a few fields once or 
twice per night. The detection rate per field and the quality of the data would be high, 
but the number of events discovered per year would be smaller, since fewer stars are being 
monitored. The total number of events discovered per year might be comparable to the 
discovery rate of the present MACHO team. 

Below we focus on modifications of the present strategies that would increase the 
detection rate of planets in wide orbits. 

9.2.3. Useful Modifications: The Monitoring Teams 

(1) Spike Analysis: Along the direction to the Bulge it seems likely that the large 
majority of the events already detected are due to lensing by ordinary stars. If a significant 
number of stars have planets in wide orbits, then short events due to lensing by these 
planets arc expected. A spike analysis of existing data from the Bulge could be productive 
in either discovering evidence of short events (in addition to those already discovered) or 
placing upper limits on the number of short-time-scale events that might have occurred. As 
more data is collected by a larger number of monitoring teams, systematic spike analyses 
should become standard. 

(2) Software Algorithms: It is important to search the data sets for repetitions, perhaps 
in conjunction with a spike analysis. 

(3) Frequent Monitoring of Some Fields: The importance of the quest for short 
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time-scale events suggests that it might be worthwhile for the teams to each choose one or 
two fields which they attempt to observe two times per night. This would allow them to 
call alerts relatively early for Jupiter-mass planets, and to have a better chance of finding 
evidence of Earth-mass planets. In addition, cooperation among the teams could greatly 
enhance the probability of finding reliable evidence of short events. At a recent workshop 
on microlensing, people associated with three of the teams discussed the advantages of 
choosing one or two fields that all of the teams would attempt to monitor each night. 
A motivating factor for such cooperation is to better understand the relative detection 
efficiencies of the teams. The quest for short events provides another important motivation. 
Indeed, if each team visited one or two Bulge fields twice per night, there would at least 
occasionally, be 6-8 times per night that those fields were checked. If, on average, 10 events 
per year were to be discovered in those fields, then, over the course of 3 years, short events 
could be discovered or ruled out at the ~ 3% level. In addition, the teams should call alerts 
for short-duration events, even if the event has apparently ceased before they can announce 
its discovery. If the short-duration event is due to lensing by a planet, calling the alert will 
allow the follow-up teams to have a better chance of detecting any subsequent repeat that 
might be due to lensing by another object in the planetary system. 

(4) Use of Pixel Techniques: During the past three years, pixel techniques have begun 
to be used for the study of microlensing in M31 (Tomaney & Crotts 1996; Crotts & 
Tomaney 1996; Crotts 1996; Ansari et al. 1997; Han 1996). It has been estimated that, 
were such techniques to be applied to the LMC and Bulge fields, the rate of event detection 
would increase by a factor of ~ 2 — 3 (Crotts 1997; Kaplan 1997). The reason for the 
increase is that observable events can occur when the baseline fiux we receive from a star 
is not bright enough for the star to appear on the templates presently used by the teams, 
if the flux is brought above the detection limits through microlensing. Such events are 
presently missed by the monitoring teams. This increase in detection efficiency, would be 
helpful, particularly in any fields singled out for frequent monitoring. Indeed, the MACHO 
team is presently engaged in applying pixel subtraction techniques to their LMC data in 
an attempt to discover events that were missed by their standard methods of detection. 
Application of these techniques to the Bulge fields would be also helpful. 

9.24. Useful Modifications: The "Follow-Up" Teams 

(1) Repeating Events: The follow-up teams need to continue to monitor all events after 
the flux has apparently fallen back to baseline. Note that this includes resonant lensing 
events, since (a) if the observed event was due to a planet, there may well be a second 
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planet, and (b) if the observed event was due to a binary lens, one or both members of the 
binary (or even the combination) may support a planetary system. 

We note that, even when the monitoring teams have identified a short-duration event 
that appears to have ceased before more frequent monitoring could begin, it is worthwhile 
for the follow-up teams to continue to monitor the flux. If the short-time-scale event was 
due to lensing by a planet, the result of frequent monitoring could be the detection of 
a repetition due to a moon revolving about the planet. If a moon exists and if we were 
fortunate in the orientation of the source track, then a repetition could occur and should be 
detected within a matter of days. A more delayed repetition could be due to later lensing 
by the central star or by another planet in the system. If they can detect events with A^ax 
as small as 1.06, the follow-up teams will be at least twice as likely to discover such events 
as are the monitoring teams. 

It should be possible for the follow up teams to discover repetitions in most cases in 
which a star is the first lens encountered. For the planetary systems we have simulated, 
roughly 1/2 — 3/4 of all other repetitions could be detected if monitoring were continued 
for 100 days. In addition, overlap events would not require a significant wait-time, since 
the second encounter would start (i.e., be associated with A > 1.06) even before the first 
encounter ended. 

(2) Isolated Events: While monitoring known ongoing events, the follow-up teams 
have many other stars in their field of view. When the total number, including those 
not individually above the detection limit (but which could be brought above the limit if 
magnified by some reasonable amount) is large enough, the follow-up teams can hope to 
identify new events. Such observations would play a unique role in the identification of new 
events, particularly the short-duration events that should be associated with wide planetary 
systems. Programs that would allow the so-called follow-up teams to take the lead in event 
detection are already underway or are planned (Sahu 1997; Gould 1997). Indeed, an ideal 
microlensing search for the purposes of the detection of planet lenses, is one in which the 
follow-up teams play the role now played by the monitoring teams to discover events, and 
continue the work by monitoring the light curves of the events they discover. 



10. Expectations 

What results are likely to be derived if the strategies sketched in the previous section 
are utilized? It is difficult to answer this question, because we know so little about how 
common planetary systems are, and about the distributions of planetary masses and orbital 
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periods. Because, however, planetary systems have begun to be discovered in the Sun's 
local neighborhood, it is beginning to seem likely that many, perhaps most stars support 
planetary systems. 

If this is so, then the preliminary results represented by the known planetary systems 
provide encouragement that microlensing will play an important role in the discovery 
of planets. The systems listed in Table 1 represent roughly 1/5 of all of the confirmed 
planetary and binary brown dwarf systems. Each would be a good candidate for detection 
via microlensing. PSR B1620-26 and the brown-dwarf system Gl 229 would produce larger 
numbers of isolated events than our solar system, but fewer repeating events. In addition, 
55 Cnc and HD 29587 are excellent candidates for lenses that would lead to resonant events. 
Q We note that, since planets in wide orbits, and even those in the zone for resonant 
lensing, are certain to be under-represented in our present census of planets, inferences 
based on Table 1 may be conservative. Thus, our present knowledge of planetary systems 
makes it seem likely that the observing teams will observe some events over the next few 
years. 

More important than the detection of any individual event, however, is our ability 
to extract information about the population of planetary systems in the regions surveyed 
by the microlensing teams. We would like to know the answers to basic questions: what 
fraction of stars have planetary systems? What are typical numbers of planets in a single 
planetary system? What are the distributions of planet masses and orbital periods? 

If the teams begin to discover planets, they will likely be able to establish the statistics 
of and distributions of properties among planetary systems in the Bulge. §10.3 is devoted 
to exploring this issue in more detail. 

10.1. Other Applications 

Before continuing with the study of planets, it is worthwhile to consider what the 
teams using the wide-planet-search strategies will learn, even if nature has been so unkind 
as to neglect to provide most stars with planets. Planet-motivated investigations will yield 
interesting fruit regardless of the size of the population of planets. 

First, they will increase the detection rates of all events, particularly short-duration 
events. Second, the frequent monitoring with good photometry will allow us to learn more 



^^47 Uma has an estimated orbital separation slightly less than 0.8 Re for D = lOkpc, x — 0.9. It may 
also be a good candidate for detection. 
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about each event detected. Common astronomical effects, such as stellar binarity, blending, 
and finite-source-size, are expected to significantly affect the shape of lensing light curves, 
introduce time dependence into the spectra, and even produce apparent repetitions. (See 
e.g, Gricst & Hu 1992; Mao & Paczyhski 1991; Mao & DiStefano 1995; DiStefano & 
Esin 1995; Kamionkowski 1995; Loeb & Sasselov 1995; Simmons, Willis & Newsam 1995; 
DiStefano & Mao 1996; Sasselov 1997; DiStefano & Perna 1997; DiStefano 1997) Studying 
the manifestations of these effects in the data sets can (1) can break the degeneracy of 
individual light curves, and (2) allow us to learn more about the populations of sources and 
lenses. Along the direction to the Bulge, we will be able to learn a good deal about the 
stellar luminosity and mass functions, and the binary fraction, as well as the distribution of 
binary properties. The information we collect can inform our design of the next generation 
of microlensing observations. For example, satellite projects have been proposed and 
preliminary calculations indicate that they are likely to be productive. (Sec, e.g., Boutreux 
& Gould 1996 and references therein.) The detailed observations that would be made 
as part of the search for planets in wide orbits would provide solid input, useful for the 
detailed planning needed for such space-based projects. 



10.2. The Relative Numbers of "Resonant" and "Wide" Events 

Without knowing more about the distributions of planet properties- which is exactly 
what we are trying to learn about through the proposed microlensing studies- it is not 
possible to make definitive predictions for the relative numbers of events due to planets 
in the zone for resonant lensing, and events due to planets in wide orbits. We do know 
enough, however, to understand the issues that determine the relative event rates. 

If we place a Jupiter-mass planet in a resonant orbit, the chance of detecting evidence 
of the planet's presence is close to 20% The present observing set-up is optimized to discover 
this type of event. If we place the same planet in a wide orbit, the probability of detecting 
an isolated short-duration event due to the planet is smaller, approximately {3n)% (where 
n is defined by Eq. 4), and the probability of detecting a repeating event in which the 
central star serves as the other lens is (6n^%)/(7ra). Thus, with the present observing 
strategy, individual examples of each type of wide-orbit event arc much less likely to be 
observed than a resonant event. Even with the present set-up, however, planets in wide 
orbits may be detected at a rate comparable to or even larger than the rate of detecting 
planets in resonant orbits, simply because on average there may be on the order of 10 times 
as many of them as there are planets in resonant orbits. 

Two factors can enhance the relative probability of detecting planets in wide orbits. 
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The first is that improved sensitivity to short-duration events and better photometric 
sensitivity can significantly increase the detection rates for all wide-orbit events, particularly 
of repeating events. The second is the influence of finite-source-size effects, which can also 
effectively increase the value of n, even for a fixed value of Amin- 

Table 4 illustrates the situation for our solar system, and for a solar system composed 
of Earth-mass planets. (See also Figure 15.) The first three lines pertain to a planetary 
system identical to our own solar system, placed in the Bulge (see also Table 3). Since, 
averaging over angles, Jupiter has a 20 — 25% chance of being viewed in the zone for 
resonant Icnsing, and since there is a ~ 20% chance of detecting a Jupiter-mass in the zone 
for resonant lensing, we have estimated that there is a ~ 5% chance of detecting evidence 
of the solar system through a resonant lensing event. For the detection criteria of set A, 
the resonant-lensing signature would be the dominant mode of detection. Changing the 
detection strategy to allow Amin — 1-06 {n — 2), allows wide-orbit events to dominate, with 
overlap and repeating events having a combined detection rate of 6.3% (7.3% if we eliminate 
the requirement of a 1— day duration for the first encounter). Thus, overlap and repeating 
events are competitive with resonant lensing events, even if we assume that we may catch 
most of the overlap events but only 1/2 of the repeating events. Isolated planet-lens events 
would be found at a rate of 1.4% and 3.1% for detection criteria B and C, respectively. 
These results indicate that the conjecture that all stars have planetary systems similar to 
our solar system, predicts the detection of resonant, repeating, and short-duration isolated 
events in roughly equal numbers. If, however, all of the planets in the model solar system 
were of Earth mass, the detection rate would fall. The fall is much more precipitous for 
resonant events, which are 5 times less likely than the combined rate of repeating and 
overlap events and 8 times less likely than isolated short-duration events. 

For the power-of-3 model, we have assumed that the probability of finding a planet 
in the zone for resonant lensing ranges from 50% to 100%. The general pattern of relative 
rates is similar to that for the solar system. For Jupiter-mass planets, resonant events are 
competitive and can even be dominant when the detection criteria of set A are used. The 
detection criteria B and C increase the detection rate of overlap and repeating events. 
Isolated short-duration events can also be important, occurring at a higher than resonant 
events. Finite-source-size effects decrease the overall detection rates, making the repeating 
and overlap events as common as resonant events, and increasing the relative importance of 
isolated planet-lens events. 

The general pattern of relative rates illustrated by these examples is likely to be 
reflected in our data sets. That is, when flnite-source-size effects are not important, the 
wide-orbit discovery channel we have studied is competitive with the resonant-event channel. 
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Table 4. Detection Rates for Planets in Resonant and Wide Orbits: 
Point and Finite-Sized Circular Sources 



Detect W 


P (2) 


Pi - Pe 


poverlap (4J 




P2 


Solar system, 


V = Gaussian, 


a = uniform: 








A 


5 


0.3 


2.1 


0.7 




B 


5 


1.4 


4.2 


1.8 




C 


5 


3.1 


4.2 


2.1 





Solar-like system w/finite source, V = Gaussian, a = uniform: 



C 03 2^4 07 1.0 

Power-of-3 system, V = Gaussian, a = uniform: 

A 10-20 13 1.9 1.5 

B 10-20 29 3.6 3.1 

C 10-20 30 3.6 3.1 
Power-of-3 system w/finite source, V = Gaussian, a = uniform: 

C 0.7-1.3 4.9 0.7 0.5 



•^^^ Descriptions of the detection conditions can be found in the text. All probabilities are given as percentages of the number 
of events in which the central star was encountered by itself, and in which the magnification reached at least A^i„ = 1.34. 

We have have used the results of Bennett & Rhie (1996) to compute the probability that a planet located there will be 
detected. They considered a mass ratio of 10~^, which, corresponds to an Earth-mass planet if the central star has a mass of 
~ O.3M0. Thus, when we use the results from Bennett $ Rhie (1996) to derive the number of resonant events expected in our 
models, we may be making a slight overestimate. To be in line with Amin for the wide-orbit calculations, we have assumed 
that the observed resonant-event perturbation should be a 6% deviation from the underlying light curve, and have interpolated 
between the results for the computation for 4% and 10% effects (Bennett $ Rhie 1996). We have also averaged uniformly 
weighted contributions for all of the cases they considered, which correspond to r ~ 1, 2, 4.3, 10, even though for the wide-orbit 
case, where we effectively also performed a linear average, we required 3 < r < 7. For the solar system, we have assumed that 
there is a 25% chance that one of the planets is in the zone for resonant lensing; for the power-of-3 model we have assumed 
that there is a 50 — 100% chance that a planet is in the zone for resonant lensing. 

Percentage of isolated (non-repeating) events (one "peak" in the light curve) with no overlap (see below) . 

Percentage of isolated (non-repeating) overlap events. In this case the two lenses were the central star and the innermost 
planet. 

Percentage of events with one repetition (two "peaks"). Overlap events are not included. 
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Repeating and overlap events tend to occur at a lower rate than either resonant or isolated 
short-duration events, but the rates of the latter two types of events are comparable. When 
finite-source-size effects are important and/or when a strategy to optimize the discovery 
of planets in wide orbits is implemented, discovery of wide-orbit planets may dominate. 
Isolated short-duration events should be the planet-lens events most frequently detected, 
and the combined rate of repeating and overlap events may be comparable to or, depending 
on the influence of finite-source size, somewhat larger than the rate of resonant- zone lensing 
events. 

10.3. The Populations of Planetary Systems and Low-Mass MACHOs 

10.3.1. Planetary Systems: General Considerations 

We found in §4 that we could hope to test simple hypotheses about specific types 
of planetary systems, even within the next few years. For example, do most stars have 
planetary systems similar to the solar system? or similar to the power-of-2 or power-of-3 
models? Of course it is most likely that planetary systems come in several different varieties. 
We can use the data to systematically extract information about lenses with a possibly 
complicated distribution of planetary system properties as follows. 

Repeating and resonant events each allow us to determine specific features of the 
planetary system that served as a lens. Each provides the value of the projected separation 
between one planet and another object (either the sun for resonant events [and most 
repeating events], or a second planet for some repeating events) in the planetary system. 
Each allows the derivation of a mass ratio. Q If more information is available, derived 
from evidence of blending or finite-source-size effects, for example, it can even be possible 
to place reasonably tight limits on the mass of the planet-lens. If we are able, therefore, to 
discover and analyze a number of repeating and/or resonant planetary-system- lens events, 
we will, by also including the influence of observational selection effects, be able to derive 
some of the characteristics of the population of planetary systems in the Galactic Bulge and 
elsewhere. 



The extraction of planet parameters, the degeneracies in those parameters, and the degeneracies between 
the planet-lens interpretation and other effects have been studied for planets in the zone for resonant lensing 
by Gaudi & Gould (1997) and Gaudi(1997). Parameter extraction should generally be more straightforward 
for planet-lenses in wide orbits. If the lens separations are very large, then the light curves are much like 
the standard point-lens light curves; if the separations are smaller, the fit given by DiStefano & Mao (1996) 
can be applied. 
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Consider a population of planetary systems that act as lenses. Let Vresi''^) be the 
probability that a planet of mass m is in the zone for resonant lensing, and Vw{Tn,a) be 
the probability that a planet of mass m is in a wide orbit, with separation a from the 
central star. As described above, discovering and measuring the rates of repeating and 
resonant events constrains the form of these probability functions. Furthermore, there is a 
consistency check, since the values of Vresip^) and P^(m, a^) must match. The integral, 

daVwiiTi, a), then predicts the rate of isolated short-duration events that should be due 
to lensing by planetary systems. 

1 0. 3. 2. Low-Mass MA CHOs 

Should the predicted rates of detectable short-duration events be larger than the 
numbers actually observed, then the analysis of planetary system lenses must be carefully 
reconsidered. We expect, however, that if there are deviations from the predictions, they 
should be because the observed rate is higher than the predicted rate. This is because 
planets in planetary systems are not the only lenses expected to produce short-time-scale 
events. For example, orbital dynamics may lead some planets to be ejected from their 
home planetary system. If even 10% of massive planets are so-ejected, then, considering 
that the microlensing observations are sensitive to the planet debris of many generations of 
stars, some of these ejected planets could be detected. In addition, it is certainly possible 
that a significant fraction of Galactic dark matter exists in the form of low-mass objects. 
It is an interesting fact that, without the information provided by repeating and resonant 
events, it will not be possible to determine whether there is a low-mass MACHO component 
in directions in which lensing by stellar systems contributes significantly to the rate of 
microlensing. Thus, when stars contribute significantly to the rate of lensing, the detection 
of repeating and resonant events is important in helping us to learn about any Galactic 
component of low-mass MAGHOs. 

11. Summary 

We have introduced the tools needed to systematically use microlensing to search for 
planets in wide orbits around distant stars. Because these searches can be conducted as 
part of the ongoing microlensing observations, we have carried out detailed simulations to 
determine what signatures should be expected if our solar system or other known or model 
systems serve as lenses. We have also studied how different detection strategies influence 
detection efficiencies. Our results are encouraging in that they clearly indicate that a 
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systematic search for planets in wide orbits is not only feasible, but that, even over the 
short term, it can yield interesting results about the population of planetary systems in our 
own and other galaxies. 

Until now, microlensing searches for planets have concentrated on searching for planets 
that might be located in the zone for resonant lensing. Resonant lensing events are 
expected and will play an important role in our quest to learn more about planets through 
microlensing. Planetary systems, however, seem likely to exhibit enough structure, in the 
form of multiple planets, moons revolving about planets, and even belts of compact debris, 
that it is important to explore all the ways microlensing can help us to study them. We 
point out that planets in wide orbits (a > 1.5Re) (1) should exist in larger numbers than 
planets in the zone for resonant lensing, (2) also can yield distinctive signatures, and (3) 
may yield larger numbers of detectable events, particularly if finite-source-size effects are 
important. Thus, previous estimates of the microlensing detection rate (see, e.g., Peale 
1997) need to be revised upward. 

The key elements of detectability are to pick up evidence of microlensing events as early 
as possible, and to be sensitive to events even if the peak magnification is smaller than 1.34. 
The detection rate of isolated planet- lens events is proportional to n, where n — Wi/RE,i, 
and Wi represents the distance of closest approach needed for reliable detection of lens i. 
The rate for repeating events scales as n'^"'"^, where k is the number of repetitions. Just as 
is the case for the detection of planets in the zone for resonant lensing, frequent monitoring 
of ongoing events is important, although hourly monitoring may not be necessary. 

The features that will improve the rate of planet detection should be considered 
when designing the next generation of microlensing observations. In the meantime, the 
MACHO results to date indicate that the teams have the necessary capabilities to detect 
wide-orbit planets. We have suggested a set of modifications in the present detection 
strategy that could, in the short-term, significantly improve detection rates for isolated and 
repeating planet-lens events. One important component is to be alert for the possibility 
that repeating events may be microlensing events, and to continue to carry out frequent 
follow-up monitoring, even after an apparently isolated event has ceased. It is true that 
the likelihood of finding a repetition may be on the order of a few percent, but the relative 
importance of the events makes the study worthwhile. After all, microlensing itself is a 
low-probability phenomenon, but has been well worth looking for. In fact, repeating events 
must be present for a number of reasons in addition to the wide-planet connection (Griest & 
Hu 1992, DiStefano & Mao, 1996). Nevertheless, almost alone among all of the light curve 
perturbations suggested so far, the search for them has remained something of a taboo. 

An important consequence of the likely presence of planets in a stellar population being 
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studied for signs of microlensing, is that isolated short events are very hkely to be present 
at a level that can be as high (compared to single stellar-lens events) as O{10%). This 
means that any signature due to low-mass MACHOs cannot be unambiguously identified 
unless the contribution due to planets can first be quantified. (There are exceptions when, 
for example, it is known that the majority of lensing events cannot be due to ordinary 
stellar systems, or when the rate of short-time-scale events is so high that the associated 
optical depth is larger than could possibly be due to planetary systems.) Fortunately, 
studies of repeating and resonant events, and any short-duration events subject to either 
finite-source-size effects or the blending of light from the central star with that of hght 
from the lensed star, should allow the contribution of planets to short-duration, apparently 
isolated events to be quantified. This contribution can then be subtracted from the total to 
derive the magnitude of any contribution due to dark matter existing outside the realm of 
ordinary stellar systems. 

Another point we have emphasized is that, even though the planetary systems 
microlensing discovers are far away, some may nevertheless be the subjects of fruitful further 
study. Indeed, when the central star is luminous, we can hope to determine its spectral 
type. In some cases, this can help to set the mass scale for the system, and can therefore 
help us to determine the mass of the planets that served as lenses. Finite-source-size 
effects can also put constraints on the lens mass. Thus, although we will not image beach 
front property on the planets discovered via microlensing, we should not give up on the 
possibility of learning more about individual planetary systems that serve as microlenses. 
We have also pointed out that it is precisely in those systems in which a planet discovered 
via microlensing is most likely to have Earth-like conditions, that the central star may be 
luminous enough to permit further study. 

Searches for planets using microlensing should be able to extend the reach of local 
planetary searches by discovering planets in distant parts of our own and other galaxies 
and by discovering even low-mass planets orbiting at low speeds. The search for planets in 
wide orbits represents a significant extension of the ongoing microlensing searches. Indeed, 
it seems likely that planets in wide orbits will provide an important, and possibly even 
the dominant mode for the detection of planetary systems via microlensing, particularly 
Earth-mass planets. 
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